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ABSTRACT 

The Tully-Fisher relation (TFR) expresses the connection between rotating 
galaxies and the dark matter haloes they inhabit, and therefore contains a wealth 
of information about galaxy formation. We construct a general framework to 
investigate whether models based on halo abundance matching are able to re¬ 
produce the observed stellar mass TFR and mass-size relation (MSR), and use 
the data to constrain galaxy formation parameters. Our model tests a range of 
plausible scenarios, differing in the response of haloes to disc formation, the rel¬ 
ative angular momentum of baryons and dark matter, the impact of selection 
effects, and the abundance matching parameters. We show that agreement with 
the observed TFR puts an upper limit on the scatter between galaxy and halo 
properties, requires weak or reversed halo contraction, and favours selection ef¬ 
fects that preferentially eliminate fast-rotating galaxies. The MSR constrains the 
ratio of the disc to halo specific angular momentum to be approximately in the 
range 0.6 —1.2. We identify and quantify two problems that models of this nature 
face. (1) They predict too large an intrinsic scatter for the MSR, and (2) they 
predict too strong an anticorrelation between the TFR and MSR residuals. We 
argue that resolving these problems requires introducing a correlation between 
stellar surface density and enclosed dark matter mass. Finally, we explore the 
expected difference between the TFRs of central and satellite galaxies, finding 
that in the favoured models this difference should be detectable in a sample of 
r\j 700 galaxies. 

Key words: galaxies: formation - galaxies: fundamental parameters - galaxies: 
haloes - galaxies: kinematics and dynamics - galaxies: spiral - dark matter. 


1 INTRODUCTION 

The complex process of galaxy formation produces sev¬ 
eral regular scaling relations, which indicate a tight con¬ 
nection between galaxies and their host dark matter 
haloes. One approach towards modelling this connection 
uses the technique of ‘(sub)halo abundance matching’ 
(AM), which sets galaxies’ stellar masses based on the 
virial mass or rotation velocity of their haloes ( [Kravtsav 
|et al.|[2004| |Conroy et al.|[2006| [Behroozi et al.||2010| 
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Guo et al. 2010| Moster et al. 2010). Further, galaxy 


sizes may be predicted from halo spin and concentra¬ 
tion under the assumption that galaxy and halo angular 
momentum are related (e.g. Fall Efstathiou]|1980 Mo 


et al. 1998). Models with these basic ingredients have 


been shown to successfully reproduce observed satellite 
fractions and galaxy two-point statistics (e.g. 
al. 2006 Reddick et al.|[^013| [Hearin et al. 


Conroy et 


2014), to 


roughly match the normalization and slope of the stel¬ 
lar mass-size relation (MSR; |Mo, Mao fc White 1998 


Kravtsov 2013), and to provide broad-brush agreement 
with the luminosity ( Trujillo-Gomez et al.|[20TT ), stellar 
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2 H. Desmond and R. Wechsler 


mass (Dutton et al.|20l0 |2011[) and baryonic (Desmond 


2012) Tully-Fisher relations (TFRs). 


In this work, we push models of this type further by 
asking whether they are capable of reproducing the de¬ 
tailed properties of the stellar mass TFR and MSR for a 
homogeneously selected sample of nearby spiral galaxies. 
Our goal is to assess whether and under which conditions 
a wide range of galaxy formation models are consistent 
with observed galaxy sizes and dynamics, and further to 
understand to what extent these observables bring addi¬ 
tional constraining power to this class of models. We ex¬ 
amine whether AM models can match the intrinsic TFR 
and MSR scatter in addition to slope and normalization, 
compare the predicted and observed correlations between 
the residuals of the two relations, and explore the ques¬ 
tion of whether satellite and central galaxies should be 
expected to lie on different TFRs. Our analysis eluci¬ 
dates and quantifies several issues which have been re¬ 
garded in some studies as fundamentally problematic for 
all A Cold Dark Matter (ACDM)-based TFR analyses 
(e.g. McGaugh||2012a|b| ). 


Most previous work on modelling the TFR has em¬ 
ployed one of three broad methodologies: semi-analytic 
modelling, AM, or hydrodynamical simulation. In semi¬ 


analytic models (e.g. Fall & Efstathiou| 1980 

iWhite & 

Frenk||I99Il [Kauffman, White & Guiderdoni| 

19931 Mo, 

Mao & White||I998| |Baughj2006| ISomerville et al.||2008| 

Guo et al. 2010 Lu et al. |20II), the relations between 


galaxy and halo properties are determined observation- 
ally or given simple parametrized forms in accordance 
with analytic calculations or the results of iV-body simu¬ 
lations. AM instead calibrates the stellar mass-halo mass 
relation directly using the observed stellar mass func¬ 
tion (SMF) and the halo mass function found in simula¬ 
tions. Assuming that the true galaxy-halo connection lies 
within the space of models accessible to this methodol¬ 
ogy, AM greatly reduces the theoretical parameter space 
and hence increases the constraining power of any set of 
observations (see also Section]^. Finally, in hydrodynam¬ 
ical galaxy simulations one attempts to account for the 
effects of baryonic physics (such as feedback from star for¬ 
mation and AGN) in a cosmological context (e.g. Abadi 


et al.|2QQ'3l|Governato et al.|20Q7|[Tissera et al.|20io|]p' 

ontek fc Steinmetz|20lH|Vogelsberger et al.|20I4| Schaye 


et al.|20f5 ). In principle, this removes the need for a priori 

assumptions concerning the relationship between galaxies 
and haloes. However, simulations of this kind are compu¬ 
tationally expensive, with the result that it is not feasible 
to model at sufficient resolution a large number of galax¬ 
ies spanning a range of values of the parameters relevant 
for the TFR. In addition to prohibiting statistical tests 
of goodness of fit, this compromises one’s ability to pin 
down the simulated relation’s intrinsic scatter. 

Our analysis expands in several ways on similar AM- 


based TFR studies in the literature (e.g. Dutton et at 
2010 Trujillo-Gomez et al. 2011). First, we consider a 


very general set of AM models, including those that have 


been shown to be consistent with galaxy clustering mea¬ 
surements (Reddick et al. 2013). This allows us to test 
the AM framework itself, not just specific AM models. 
Secondly, we include in our framework a prescription for 
taking account of the impact of selection criteria in the 
observational studies, which to first order eliminate non¬ 
rotating early-type galaxies. Thus we intend our theoret¬ 
ical TFR and MSR populations to contain only galax¬ 
ies that would pass the cuts associated with the data. 
Thirdly, in addition to considering the slope, intercept, 
and scatter of the TFR, we also investigate the shape of 
the stellar MSR and the correlation between the residu¬ 
als of these two relations. This information provides en¬ 
hanced constraining power over the TFR alone. Finally, 
we evaluate our models in a statistically rigorous way by 
means of a full likelihood formalism. 

The structure of this paper is as follows. Section 
describes the observational data with which we compare 
our model, and Section [^documents the iV-body simula¬ 
tions from which halo properties are drawn. Section [^de¬ 
tails our methodology. In Sectionwe present our results 
for the TFR, MSR, correlation between the residuals of 
the TFR and MSR, and environment dependence of the 
TFR. In Section we discuss the general significance and 
broader ramifications of our results. Section proposes 
avenues for future research, and Section [^concludes. Ap¬ 
pendix provides more specific information on the ef¬ 
fect of the individual model parameters on the TFR and 
MSR. In Appendix [B] we discuss in detail the assump¬ 
tions in our model and the systematic errors that may 
result from their failure, and document the tests we have 
performed to assess them quantitatively. We provide ad¬ 
ditional technical implementation details in Appendix [C] 


2 OBSERVATIONAL DATA 


The observational TFR to which we compare our model 
is that of Pizagno et al. (2007 hereafter POT). This is 
a sample of 162 galaxies with Ha rotation curves ob¬ 
tained via long-slit spectroscopy from the Galar Alto and 
MDM observatories, and luminosities measured by the 
Sloan Digital Sky Survey (SDSS; York et al.|2000 ). The r- 
band magnitude distribution is roughly uniform between 
-22 and -18.5, corresponding to stellar masses roughly in 
the range 10^ — Mq. The sample was selected from 
SDSS DR2 according to four criteria: 


(i) cz > 5000 km s“^. This ensures that peculiar veloc¬ 
ities are small relative to the Hubble flow. 

(ii) cz < 9000 kms“^ for —18.5 > Mr > —20, 
IlOOOkms”^ for —20 > Mr > —21, and 15, 000 km s“^ 
for Mr < — 21. This eliminates distant galaxies with more 
uncertain velocity measurements. 

(hi) An isophotal axis ratio cut hja < 0.6, where h is 
the apparent minor axis of the galaxy and a the appar¬ 
ent major axis, as measured in the r band by the SDSS 
pipeline. Focusing on edge-on galaxies in this way reduces 


© 0000 RAS, MNRAS 000 , 000-000 
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uncertainties in inclination measurements and increases 
the likelihood that a given galaxy is a rotation-dominated 
spiral. 

(iv) The rotation curve is well fit by an arctan func¬ 
tion, as determined by visual inspection. This eliminates 
morphologically unusual galaxies and those without sig¬ 
nificant rotation. 


As we describe further in Section we take the axis 
ratio cut (which we can systematically investigate in the 
SDSS) to be the selection criterion that most strongly 
influences the position of the accepted galaxies on the 
TFR plane, and assume that the arctan cut eliminates 
a random subsample of these galaxies (i.e. whether or 
not a galaxy passes the latter cut is independent of the 
properties of its host halo that influence its position on 
the TFR). Further discussion of this assumption can be 
found in Appendix [B] 

The POT study itself uses luminosity; the conversion 
to stellar mass requires a mass-to-light ratio for each 
galaxy, in addition to an extrapolation of the observed 
surface brightness profile to find the total luminosity. Sev¬ 
eral different prescriptions for these mappings exist in the 
literature, sometimes based on quite different assump¬ 
tions and underlying stellar population models (e.g. [Bell 


fc de Jong|200H|Kauffmann et al.|2003|[McGaugh|2005b| 

Bernard! et al.|2013| ) , making this step a significant source 
of systematic uncertainty in most studies of the stellar 
mass TFR. In this work we take the stellar masses of 
the POT galaxies from the NASA Sloan Atlas (NSA). 
This data set uses improved sky subtraction optimized 
for low-redshift galaxies, and is in good agreement with 


recent estimates of the SMF from Bernard! et al. (2013) 


and Kravtsov et al. 


(2014) 


The advantages of the POT sample for the purposes 
of our study are threefold. First, the selection criteria are 
relatively loose and simple to model theoretically, which 
is important because selection effects are typically one of 
the largest (and least discussed) sources of uncertainty 
in comparison of data with models that make predictions 
for the full galaxy population. As described in Section 
we model the impact of selection by postulating and then 
marginalizing over a correlation between a galaxy’s axis 
ratio and rotation velocity. Secondly, the fact that the 
galaxies were selected for spectroscopic follow-up from a 
well-understood parent population (those in the SDSS) 
allows us to achieve consistency between the simulated 
and observed galaxies using correlations from the SDSS. 
Finally, the POT study was done with the aim of theory 


^ We remove five POT galaxies (J13213.08+002032.T, 
J012223.T8-005230.T, J124545.20+535T02.0, 

J1T0310.4T+65341T.6, J20530T.50-00240T.0) whose 

masses are not recorded in the NSA, and another 
(J215652.T0+12185T.5) with stellar mass significantly 
lower than the rest. This leaves 156 galaxies in our TFR 
sample. 
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comparison explicitly in mind. Thus, for example, they 
use a velocity measure that is easy to model robustly 
(the velocity at the radius enclosing 80 per cent of the 
i-band light, Vso), and expend considerable effort on pin¬ 
ning down the relation’s intrinsic scatter. 

For completeness, we note that a study similar to 
POT was more recently performed by Reyes et al. (2011). 
We choose to use the former because its selection criteria 
are simpler and hence easier to model within our frame¬ 
work. The two data sets are sufficiently similar that our 
conclusions would unlikely change, were the [Reyes et al.| 
(2011) data used instead. Other sizeable TFR samples in 
the literature - which do not possess the homogeneity or 
simple selection required for our detailed modelling - in¬ 
clude those of Courteau et al. ( |200T| , [McGaugh (2012b), 
and jMocz et al. (2012). 

Although the baryonic TFR (which includes cold gas 
mass) appears to have smaller scatter and curvature than 
the stellar mass TFR ( McGaugh et al.||2000 ), and hence 
to be more fundamental, practical considerations force us 
to restrict our attention to the latter. Homogeneous data 
sets that both contain cold gas mass measurements and 
are derived from a thoroughly characterized parent pop¬ 
ulation do not yet exist. Stellar masses can be estimated 
from the information acquired by large-scale galaxy sur¬ 
veys such as the SDSS, whilst cold gas mass measure¬ 
ments (requiring Hi observations) have been done for rel¬ 
atively few systems not systematically selected according 
to a fixed set of criteria. A further consequence of this 
is that the SMF, required for AM, is known with signif¬ 
icantly better accuracy than the corresponding baryon 
mass function. An analysis at our level of detail is not 
yet feasible for the baryonic TFR, but will be worth per¬ 
forming when it becomes so. 


3 SIMULATION DATA 


The choice of A-body simulation with which to create 
our theoretical galaxy population is driven by two com¬ 
peting criteria. On the one hand, we require the sim¬ 
ulation to contain a statistically representative sample 
of haloes hosting galaxies with stellar masses at the 
upper end of the range covered by the POT data set 
[log(M*/M 0 ) 11 . 2 ]. On the other, we require haloes 

at the low mass end [log(M*/Mo) 9.0] to be suffi¬ 
ciently well resolved for their concentrations to be reli¬ 
able. Diemer fc Kravtsovj ( |2015| propose three criteria 
for a halo to be well resolved. 


(i) At least 1000 particles within R 200 C, the halo radius 
within which the mean density is 200 times the critical 
density of the Universe. 

(ii) At least 200 particles within rg, the scalelength of 
the halo. 

(hi) Ts at least 6 times the force softening length, e, 
used in the simulation. 
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Many log(M*/M0) 9.0 haloes (using a fiducial 
AM prescription) fail at least one of these cuts in a 
(250 Mpc h~^)^ simulation with 2048^ particles and e ~ 
1 h~^ kpc. We find that this failure is manifest in a > 
10 per cent discrepancy between median concentrations 
from the simulation and those predicted by the model 
of Diemer & Kravtsov (2015), and > 10 per cent differ¬ 
ences between alternative methods of determining con¬ 
centration. Thus we must use a higher resolution simula¬ 
tion to accurately model low-mass galaxies. To this end, 
we use the c-125 simulation, a (125 Mpc box with 

2048^ particles and comoving softening length 0.5 h~^ 
kpc, simulated using l-GADGET (based on gadget- 
Springel|2005 ). This simulation { 


2; Springel et al. 


2001 


sumes a flat ACDM cosmology with h = 0.7, = 0.29, 

Qh = 0.047, (78 = 0.82 and Us = 0.96. The initial condi¬ 
tions were generated by 2lptig ( Crocce et al.|20Q6 ) at 
= 199, and the power spectrum was generated by GAM^ 
The c-125 box, however, contains few high-mass 
haloes due to its relatively small volume. For exam¬ 
ple, there are only ~ 220 haloes in the simulation that 
would be assigned galaxies with 11.15 < log(M*/M0) < 
11.2 by a standard abundance match - too few to jus¬ 
tify statistical inference concerning the high-mass end 
of the TFR. Modelling this regime accurately therefore 
requires a larger box. We use a larger simulation run 
with the same code and cosmological parameters for this 
purpose. The c-250 box contains 2048^ particles in a 
volume of (250 Mpc yielding a particle mass of 

I. 46 X 10^Mq and has a comoving softening length 
of 1 h~^ kpc. This box contains ~ 1800 haloes with 

II. 15 < log(M*/M0) < 11.2. Both boxes were provided 
by M. Becker (Becker et ah, in preparation). We identify 
haloes using the ROGKSTAR halo finder ( [Behroozi, Wech- 


sler & Wu 2013a), and generate merger trees using the 
GONSISTENT TREES algorithm ( Behroozi et al.||2013b ). 

To create a final galaxy-halo catalogue with which to 
model the TFR, we proceed as follows. We first perform 
AM separately on the haloes generated by the c-125 and 
c-250 simulations (see also Section]^. We then splice to¬ 
gether the log(M*/M0) < 10 part of the c-125 catalogue 
and the log(M^/M0) > 10 part of the c-250 catalogue. 
We have verified that the median halo concentrations, 
virial masses, and spins of the two simulations differ by 
< 5 per cent for 9 < log(M*/M0) < 11.2, so the discon¬ 
tinuity at log(M*/M0) = 10 is not significant. Our final 
catalogue contains at least 1800 objects in all stellar mass 
bins, and all haloes are fully resolved. 

In order to depend as little as possible on the as¬ 
sumption that haloes have an NFW density profile, we 
use Mvir values determined directly by the halo finder 
(rather than calculating them from rvir), and determine 
concentration from the ROGKSTAR output Vs,kiypin’ (de¬ 
rived from the ratio of r’max to r'vir) rather than rg. 


^ http://camb.info/ 


4 METHOD 

We begin with an executive summary of our method for 
comparing mock and observed TFRs and MSRs. Each 
step is expanded upon in the remainder of this section. 

(i) Use AM to associate a stellar mass with each halo, 
using one of six halo properties as the AM parameter and 
allowing for the possibility of non-zero scatter. 

(ii) Assume the stellar mass is distributed in a thin 
exponential disc and find the scalelength by requiring the 
specific angular momentum, J, of the disc to equal some 
fraction of the specific angular momentum of the host 
halo. 

(hi) Model the effect of disc formation on the density 
profile of the halo, allowing for standard adiabatic con¬ 
traction, an expansion of the same magnitude using the 
same functional form, or anything in between. 

(iv) Add a bulge to each model galaxy by bootstrap 
resampling from the bulge-to-disc ratios of the P07 galax¬ 
ies, in five bins of stellar mass. 

(v) Model the impact of the selection effects in the 
P07 data set by imposing a correlation between apparent 
axis ratio and Vso at fixed stellar mass. The strength 
of this correlation is controlled by a free parameter that 
interpolates between no correlation and a near-monotonic 
relationship. 

(vi) Generate 200 mock data sets containing as many 
galaxies as are in the P07 sample, and with the same stel¬ 
lar masses and velocity uncertainties, but with velocities 
chosen randomly from the theoretical parent population. 
Fit each data set with a power-law TFR to determine 
the best-fitting slope, intercept, and scatter values, and 
compare the distributions thereby obtained to the values 
for the real data in order to assess goodness of fit via a 
likelihood. 

(vii) Use an MCMC (Markov Chain Monte Carlo) al¬ 
gorithm to explore the parameter space systematically 
and obtain posterior probability distributions for the 
model parameters. 

(viii) Repeat the above set of steps replacing the TFR 
by the MSR or the correlation of the residuals of the two 
relations. 


Our model is built on the technique of AM ( [Kravtsov et 
al.||2004[ [Conroy, Wechsler Kr avtsov|| 2006| [Behroozi, 


Conroy & Wechsler 2010| Guo et al.[[2010| [Moster et 


al. 


2010 


Reddick et al. 2013). This enables the asso¬ 


ciation of a galaxy with each dark matter halo in an 
Wbody simulation by hypothesizing a one-to-one rela¬ 
tion between stellar mass and some halo property. One 
begins by rank-ordering the haloes in the simulation at 
= 0 by that property, and then uses a SMF (in our 
case that of Bernardi et al.[[2013 ) to ascertain the abun¬ 
dances of galaxies as a function of stellar mass that one 
would expect to see in the volume of the simulation. 
One then rank-orders these galaxies by stellar mass and 
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enforces the galaxy-halo connection by associating the 
galaxy with nth largest stellar mass to the halo with 
the nth largest value of the matching property. AM has 
been shown to produce galaxy catalogues that are con¬ 
sistent with a large number of observed statistics, includ¬ 
ing galaxy clustering, galaxy kinematics, galaxy-galaxy 
leasing, galaxy group catalogues, and galaxy void statis- 


tions and galaxy two-point correlation functions. Similar 
results have been found by other authors using satellite 
velocity dispersions (e.g. More et al. 2009). We investigate 
a range of AM models rather than restricting ourselves 


to those deemed acceptable by Reddick et al. (2013) in 


order to explore in full generality the aptitude of the AM 
ansatz for reproducing the observed TFR and MSR, and 


tics (Conroy et al.|2006 

Vale V Ostriker|2006 Moster et 

to bring out important physical differences between the 

al.||2010| Reddick et al. 

2013). It has also been shown to 

models in terms of the internal galaxy dynamics they im- 


yield a stellar mass-halo mass relation similar to that ob¬ 
tained directly from studies of weak leasing and satellite 
kinematics ( Behroozi et al.|2010 ). 


The simplest halo property to use for AM is the cur¬ 
rent virial mass of the halo, Mnow Much effort has re¬ 
cently been expended, however, on generalizing the pro¬ 
cedure by considering alternative halo properties. For ex¬ 
ample, it might be expected that the process of galaxy 
formation would be more closely correlated with the 
depth of the potential well in which the galaxy is located, 
for which a better proxy than Mnow is the present-day 
maximum rotation velocity of the halo, Mnax- However, 
galaxies did not form at z = 0; the stellar mass that a 
given halo acquires may in fact be better traced by the 
peak value of the halo’s maximum rotation velocity over 
its merger history, Rpeak- Satellite galaxies, on the other 
hand, likely form when the subhaloes that harbour them 
first cross the virial radius of the host halo; thus another 
popular choice (Wcc) is to use the maximum rotation 
velocity at the time of accretion for subhaloes, and the 
current maximum rotation velocity for distinct haloes. 
Along with the analogues of Vacc and Rpeak obtained by 
replacing the halo’s velocity by its mass (Mace, Mpeak), 
these are the models that we investigate in parallel in this 
work0 

It is also possible to loosen the galaxy-halo connec¬ 
tion within the AM framework. This is done with a pa¬ 
rameter known as ‘AM scatter’, measured in dex, which 
sets the standard deviation of a normal distribution. Af¬ 
ter completing AM with zero scatter, the stellar mass 
of each galaxy is perturbed by an amount drawn from 
this distribution, and then the abundance of galaxies as 
a function of stellar mass is modified to restore consis¬ 


tency with the true SMF (see Behroozi et al. 2010 for 


details). The larger the AM scatter, the larger the scat¬ 
ter in the relationship between stellar mass and halo mass 
or velocity. We take the AM scatter to be a free parame¬ 
ter with uniform prior in the range 0-0.5 dex, and assume 
for simplicity that it is mass independent. 


The analysis of Reddick et al. (2013) suggests that 


only by matching to Rpeak with a scatter of 0.2 ± 0.03 
dex can one simultaneously fit observed satellite frac- 


^ We have also tried abundance matching to the maximum 
rotation velocity of a halo at the time it attains its maximum 
mass, Mnpeak- This generates a very similar galaxy-halo con¬ 
nection to the Fpeak model, and hence an almost identical 
TFR. 


et al. (2013) is optimal for distinguishing the models on 


the basis of their assembly bias and satellite fractions, the 
TFR is most suited to bringing out differences in the stel¬ 
lar mass-concentration relations. We use the Rpeak model 
most often for illustrating our results in recognition of the 
fact that it is currently the most favoured observationally, 
but highlight where alternative proxies differ. 

After setting the stellar mass of the galaxy within 
each halo in this way, we model its structure. We assume 
each galaxy to be a thin exponential disc, a good approx¬ 
imation for the late-type spirals under investigation. We 
set the disc’s scalelength by the requirement that its spe¬ 
cific angular momentum be some fixed fraction of that 
of its host halo (as determined by the halo spin param¬ 


eter, A), a common assumption since the work of Mo et 
and still popular today. The ratio of the disc 
and halo specific angular momentum is the third free pa¬ 
rameter in the model: j = We assume this to be 

a constant for all haloeslj Since a priori expectations for 
j from hydrodynamical simulations range from j 1 


(e.g. Navarro & White||1994| |Navarro & Steinmetz 

199T 

Scannapieco et al.|2012) to j > 1 (e.g.|Kimm et al. 

2011 

Stewart et al.|2013 

Danovich et al.|2015), we begin with 


a uniform prior over the range (0,5]. We will find in Sec¬ 
tion that j is constrained by the MSR to be < 1.4 at 
3(7 for all AM proxies, while it is poorly constrained by 
the TFR. To prevent a broad j posterior from degrad¬ 
ing the constraints on the other model parameters with 
which it is degenerate, we restrict j to the range (0,1.4] 
for our primary TFR analysis. We discuss the effect of 
the j prior further in Section |5.1| 

We next add a bulge to our model galaxies to further 
augment consistency with the POT data set. To avoid as¬ 
sumptions about the distribution of bulge-to-disc ratios 
in our parent population, and their correlation with other 
galaxy properties, we do bootstrap sampling of their val¬ 
ues for the galaxies in the POT sample. We partition 
the POT galaxies into five bins of log(M^/Mo) (9-9.5, 
9.5-10, 10-10.5, 10.5-11, 11-11.5), and then select for 


^ It is possible that j is mass dependent. However, we show 
below that no mass dependence is required in order to fit the 
slope and intercept of the TFR or MSR. A mass-dependence 
term included in the model would be poorly constrained and 
consistent with zero, and would be unlikely to alter any of our 
conclusions. 
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each model galaxy a bulge-to-disc ratio from the corre¬ 
sponding bin|^We assume that the bulge-to-disc z-band 
flux fractions measured by POT equal the corresponding 
mass fractions, and defer to Appendix [B] a discussion of 
the (small) magnitude of the systematic error that may 
thereby be induced. The bulge and disc masses sum to 
the galaxy stellar mass assigned by AM. 

The disc scalelength, Rd, that gives the galaxy the 
correct angular momentum is determined by the density 
prohle of the halo. Although we assume the pure dark 


matter haloes to have NFW density prohles (Navarro 


Frenk & White 1996), we allow for the possibility that 


these are modihed by disc formation. In particular, the 
collapse of the baryons into a relatively small galaxy 
at the centre is expected to draw in the surrounding 
dark matter in a process known as ‘adiabatic contrac¬ 


tion’ ( Blumenthal et aL]|1986 Gnedin et al.|20(M 2011). 
Our fiducial model of adiabatic contraction in this study 


is that of Gnedin et al. (2011), which builds upon earlier 


work by accounting for non-circular dark matter orbits. 
In this model 


Mf(rf) rf = Mi(n) n, (1) 

where subscripts i and f denote before and after disc for¬ 
mation, respectively. 


Mi(r) = Mi,halo (r) (2) 

is the initial mass enclosed within radius r (with the halo 
described by an NFW prohle). 


Mf(r) = Mdisc(^) + Mbulge + (1 — M*/Mvir) Mf,halo(^) 

(3) 

is the hnal enclosed mass, and 


r = Aq To {r/roY 


(4) 


where Aq = 1.6, w = 0.8, and ro = 0.03 i^vir (see Gnedin| 
et al.||2011 eq. 4). For a given halo mass, concentration, 
and spin, and galaxy disc and bulge mass, these equa¬ 
tions may be solved iteratively with the requirement that 
^disc = j e/haio to determine the hnal halo mass within 
any given radius, and the disc scalelength. This gives 
the ‘standard’ adiabatic contraction solution. However, 
it is commonly argued (largely motivated by a desire to 
solve various small-scale problems with AGDM such as 
the ‘cusp/core’ - de Blok||2010|- and ‘too b ig to fail’ - 
Boylan-Kolchin, Bullock fc Kaplinghat|20TT - problems) 
that processes such as baryon ejection due to stellar winds 
and energy injection by supernovae can in fact generate 


^ This procedure coarsely incorporates the correlation found 
in the data between bulge mass fraction and stellar mass. Cor¬ 
relations between bulge mass fraction and other properties of 
the POT galaxies are found to be negligible. 


an expansion of the halo that erases the impact of the ear¬ 
lier stage of contraction. We remain agnostic about the 
detailed effect that baryonic feedback has on halo density 
prohles and choose to employ a simple generalization of 
the adiabatic contraction formalism created by |Dutton| 


et al. (200T). In particular, we dehne 


F(ri) = rf/n, (5) 

and then take the true hnal radius enclosing a given 
amount of mass to be given by 


^/,true — F Ti. 


( 6 ) 


The free parameter d interpolates between standard adi¬ 
abatic contraction (i^ = 1), no effect of baryonic collapse 
on the dark matter halo (z/ = 0), and halo expansion by 
the same factor as the standard case (y — — l)[jDifferent 
values of v result in different disc scalelengths and halo 
density prohles for given input parameters. 

We now have all the information needed to calculate 
the full rotation curve of all model galaxies, and hence 
locate them on the TFR. To do so, we record the velocity 
at the radius which encloses 80 per cent of the Tband 
light (converting from stellar disc radius to z-band hux 
radius using a correction factor of 1.13; [Dutton et al.| 
20QT| ). Under the assumption of a constant stellar mass- 
to-light ratio over the bulge and disc (discussed further 
in Appendix H- this velocity is the model analogue of 
POT’S Uso. 


Finally, we model the potential impact of the selec¬ 
tion criteria employed by POT. Since we cannot deter¬ 
mine which of our model galaxies would have properties 
allowing them to pass the POT selection cuts (that is to 
say, we do not know how to map selection criteria on 
to halo properties), we attempt merely to bracket their 
effect on the TFR. Thus we assume that there exists a 
correlation between apparent axis ratio and the Vso value 
of a galaxy, with larger axis ratio implying larger VsoQ 
There is weak evidence for such a correlation - galax¬ 
ies with larger apparent axis ratios are more likely to be 
early type, which have been found to lie preferentially 
in more massive or concentrated haloes at hxed galaxy 
mass (e.g. Rodrfguez-Puebla et al.pQll Wojtak & Ma- 


mon||20I3 ). However, we do not motivate this model on 


® There is no a priori theoretical reason why u must lie in the 
range [-1,1]. However, we show below that the TFR constrains 
u to lie comfortably within these bounds, and hence a larger 
range need not be considered. 

^ An alternative would be to use the age matching model 
of |Hearin Sz Watson] ( [^13| ) to assign a colour to each model 
galaxy, and then assume a correlation between apparent axis 
ratio and colour. We have found that the correlations between 
halo age and the dynamical properties which affect Vso are 
too small for selection effects to have a significant impact on 
the TFR, under this assumption. 
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the strength of this evidence - we wish only to assess the 
extent to which selection effects could shape the TFR, 
and provide a framework for taking these into account 
when comparing theory with observation. Since galaxies 
are typically predicted by AM to rotate too fast for their 


stellar mass at log(M*/M0) > 10 (e.g. Trujillo-Gomez et 
al.|2011 Desmond|2012| , removing the galaxies with high 


Vso in this way is expected to improve agreement with the 
observations. Our formalism allows us to determine how 
significant such selection effects must be. 

The details of the implementation of this correlation 
in our model are as follows. We begin by determining the 
fraction of galaxies in the local Universe, as a function of 
stellar mass, which pass the POT axis ratio cut {b/a < 0.6; 
see Section]^. This is done using the SDSS Value Added 
Galaxy Gatalog (VAGG; Blanton et al. 2005), and de¬ 
termines the fraction of galaxies that should be included 
in the theoretical parent population under construction. 
The acceptance fraction falls from around 40 per cent at 


Next, we establish at each stellar mass a probability den¬ 
sity function on the set of Vso values output by our model 
at that mass, which by construction describes the likeli¬ 
hood of a galaxy with that Vso value having b/a < 0.6, 
and thus the likelihood that it should be included in our 
catalogue. This function is made linear, and the ratio 
of the probability density at the smallest Vso value pro¬ 
duced by the model at that mass to the largest is taken 
to be another free parameter in the model, with a value 
in the range [0.5,oo)|^To compactify this semi-inhnite 
parameter range into a finite one, we perform the re- 
parametrization y = arctan(x — 0.5) where x is the origi¬ 
nal parameter and y the new one. Thus y (which we refer 
to hereafter as the ‘selection factor’) lies in the range 
[0,7r/2). 

In equations, we calculate the probability of accept¬ 
ing a model galaxy with Vso = V (that is, the probability 
of assuming it has b/a < 0.6) as 


Paccept(U)=A-FBU, 

where A and B are constants set by 


(7) 



Taccept (Vso) dVsO — 1 


and 


P(Umin) ^ 
P(Knax) 

( 8 ) 


Finally, we select the number of galaxies at each stel¬ 
lar mass corresponding to the fraction determined from 


® The lower bound of 0.5 is arbitrary, and simply limits the 
possibility for strong preferential selection of high Vso galax¬ 
ies at fixed stellar mass. Such selection would conflict with the 
findings of Rodrfguez-Puebla et al. ]( |2011|> an d |Wojtak fc MaT] 
|mon| ( [20T3 1, and will be shown in Section |5.1| to be disfavoured 
by the TFR. 


the VAGG, according to this probability density func¬ 
tion. The larger the value of the selection factor, the more 
likely galaxies are to be selected from the low end of the 
Vso function at fixed stellar mass, and hence the stronger 
the preference in the model for TFR selection effects to 
eliminate high-velocity galaxies. This completes the cre¬ 
ation of a full theoretical parent population of galaxies, 
from which, absent systematic errors, the POT sample 
should have been randomly drawn. Each set of values 
for the five free parameters in our model (summarized in 
Table generates a different parent population. 

It now remains to assess the extent of agreement be¬ 
tween the theory and POT data. To do this, we generate 
a large number (200) of mock data sets from the theo¬ 
retical parent population, each one containing the same 
number of galaxies as the POT sample, and with the same 
masses and velocity uncertainties, but with velocities ran¬ 
domly selected from the model galaxies. To simulate the 
effect of measurement error, we perturb each model veloc¬ 
ity by a random number drawn from a normal distribu¬ 
tion with mean zero and standard deviation equal to the 
measurement uncertainty. We then fit a power law to the 
log(V8o/km s“^)-log(M* /Mq) relation of each mock data 
set (using a Gaussian likelihood model) and record the 
best-fitting values of the slope, intercept and scatter, for 
comparison with the analogous quantities derived from 
the observational data. We pivot the fit on the median 
log(M*/M0) of the sample, 10.4T, to eliminate the de¬ 
generacy between the slope and intercept, and take stel¬ 
lar mass to be the independent variable throughout. To 
condense this information into a single measure of good¬ 
ness of fit, we locate each mock data set in the space of 
{slope, intercept, scatter} values and fit a 3D Gaussian to 
the 200 points thereby obtained. This describes the the¬ 
oretical TFR probability density function in this space. 
The value of this function evaluated at the {slope, in¬ 
tercept, scatter} coordinate of the POT data (0.264, 2.24, 
0.0T80) is taken to be the likelihood of the model un¬ 
der consideration: the further the POT coordinate from 
the centroid of the theoretical distribution, the worse the 
agreement between theory and observation. A goodness 
of fit measure using mock data sets with characteriztics 
identical to the real data treats simulated galaxies in the 
same way as observed galaxies, and is therefore a more 
sound basis for statistical inference than a measure based 
on the entire model population. 

Thus we obtain the likelihood of the observational 
data given a particular set of model parameter values. We 
now use Bayes’ theorem to convert this into a probability 
of the parameters given the data (using uniform priors 
over the allowed ranges): 


Pip\d) oc 





(9) 


where p = {AM parameter, AM scatter, j, v, selection factor} 
is the vector of model parameters, d = {M*,V8o} is 
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Description 

Allowed values 

AM parameter 
AM scatter 

3 

u 

Selection factor 

Halo property to which galaxy stellar masses are matched 
Universal Gaussian scatter in stellar mass at fixed halo property 
Ratio of specific angular momentum of disc and halo 

Controls degree of halo expansion (-) or contraction (+) 

Governs impact of POT selection effects on theoretical TFR 

^^eak) ^^aX) -^peakj -^acci -^now 

[0, 0.5] dex 

MSR: (0, 5] TFR: (0, 1.4] 

[-C 1] 

[0, 7r/2) 


Table 1. Table of model parameters and their allowed values. 


the POT data vector, C is the 3x3 covariance ma¬ 
trix of {slope, intercept, scatter} values of fits to 
the mock data sets, a? is a 1 x 3 vector of the mean 
{slope, intercept, scatter} values of the mock data, 
and Xd = (0.264,2.24,0.0780). Finally, we calculate the 
posterior distributions of the model parameters with 
th e MCMC algorit hm. We employ the sampler created 
by Lu et al. (2014) fusing 32 walkers and at least 4000 
steps after burn-in, and check convergence by eye. Since 
the AM parameter is a discrete variable, we perform the 
MCMC analysis separately for each of its six possible 
values. 

To further assess the agreement of our model with 
the POT data, and to further constrain our model param¬ 
eters, we examine the relation between log(M*/M0) and 
log(i?d/kpc) (MSR) using the same formalism. Since the 
disc scalelengths were measured by POT (and were calcu¬ 
lated in our model en route to Vso), this simply requires 
replacing Vso by Rd in the final step and then repeating 
the MCMC. The best-fitting {slope, intercept, scatter} of 
the POT MSR is {0.233, 0.486, 0.170}. We also investigate 
the cross-correlation of the TFR and MSR by generating 
combined radius and velocity mock data, and establish 
goodness of ht using a likelihood based on the strength 
of the correlation of the radius and velocity residuals. 
Further discussion of these analyses and the results they 
yield may be found in Sections |5.2| and |5.3| 

More technical details of our method (in particular 
to make each likelihood evaluation quick enough for an 
MCMC analysis to be feasible) are given in Appendix [C| 


5 RESULTS 

5.1 Tully—Fisher relation 

Table [2] lists the Icr and 2cr constraints on our model 
parameters after comparison to the POT TFR, and Ta- 
blej^gives the maximum-likelihood parameter values and 
their associated goodness of ht. Fig. shows the poste¬ 
rior probability distributions obtained for the parame¬ 
ters, and Fig. compares the distribution of slope, in¬ 
tercept and scatter values obtained from our mock data 
sets using the maximum-likelihood parameter values to 


® https://github.com/ylu2010/PSEtoolbox 


those of the data. In Appendix [A| (Fig. Al), we show 
and explain the effect of each of the model parameters 
on the predicted TFR. The most signihcant hndings are 
the following. 

• In ah cases, reasonably low AM scatter (e.g. < 0.36 
dex at 2cr for the Upeak modef^ is required in order to 
limit the scatter in the model TFRs to the small value 
observed in the data. This constraint is currently weaker 
than that afforded by clustering and satellite fraction 
measurements (e.g. More et al. 2011| Reddick et al. 2013), 
but can be strengthened when larger homogeneous TFR 
data sets become available. In addition, the TFR has the 
potential to constrain the scatter at lower mass (see also 
Section 


6 . 1 ). 


• The Upeak model requires at the 2cr level that galaxies 
passing the POT selection cuts have systematically low 
rotation velocities for their stellar mass (see hnal column 
of Table and bottom-right panel of Fig. 1(a)). This 
further reduces the predicted scatter and also improves 
agreement with the data by shifting the TFR towards 
lower Vso at hxed stellar mass. All other AM models are 
consistent at 2cr with random selection [i.e. independent 
of Vso, corresponding to selection factor = arctan(0.5) = 
0.464]. It is important to note, however, that models with 
weak selection require high j to be consistent with the 
relatively low normalization of the TFR, and hence that 
lowering the upper limit of the prior on j would tighten 
the selection requirements. For example, if j is capped at 
1 instead of 1.4, only the Mnow and Mace models would 
be consistent with random selection at the 2cr level. 

• For ah AM models, standard adiabatic contraction 
would make galaxies rotate too fast at hxed stellar mass 
and is ruled out at the > 4a level. For mass-based AM 
models, agreement with the data is maximized by switch¬ 
ing contraction off altogether. In contrast, for models that 
match to halo velocity, expansion of the halo in response 
to galaxy formation is in fact required at the 2cr level. 


Technically, each AM scatter value presented here is the 
quadrature sum of the uncertainty in the input SMF and the 
intrinsic scatter associated with the AM procedure itself. So, 
for example, if the stellar masses used to construct the SMF 
have an uncertainty of 0.15 dex, then our Fpeak constraint 
implies that the pure AM scatter must be < 0.33 dex at 2a 
or < 0.10 dex at Icr. See Appendix [B] for further discussion of 
stellar mass uncertainties. 
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(6) (f) ^now 

© 0000 RAS, MNRAS 000, 000-000 

Figure 1. Posterior probability distributions of model parameters obtained from the TFR (diagonal panels), and correlations 
among the parameters (off-diagonal panels). Contours indicate the 68 and 95 per cent confidence levels. Each subfigure shows the 
result for a different AM model; note the change of range on the u axis between the M and V models. Light blue lines show the 
case of no effect of disc formation on halo density profile {u = 0) and random selection [selection factor = arctan(0.5) = 0.464]. 
In several places in the paper we compare the Epeak (fop left) and Mnow (bottom right) models, which span the range of possible 
effects of AM proxy on the model TFR. 
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(However, we show in Appendix that a small amount 
of contraction is permitted if Qm and as take their lowest 
values allowed by the Planck uncertainties.) We compare 
the v posterior probability distributions for the Hpeak and 
Mnow models in Fig. 

• Low j values (j ^0.7, depending on the model) are 
excluded. Lowering j reduces the angular momentum of 
all galaxies, tending to decrease their size and increase 
their rotation velocity. Our main analysis uses the prior 
0 < j < 1.4, based on constraints from the MSR dis¬ 
cussed in the following section. If this prior is relaxed, 
the upper bound placed on j by the TFR ranges from 
~ 2.5 for the Mnow model to ~ 5 for the Rpeak model. In 
this case, high u values {u > 1) are allowed, and no lower 
limit may be placed on the selection factor for any AM 
model. 

• The trends in the constraints derived using the dif¬ 
ferent AM models may be understood by means of the 
trend in average halo concentration as a function of stellar 
mass. The models in Tableare listed in decreasing order 
of average halo concentration over the stellar mass range 
of interest (with the largest difference being between Mnax 
and Mpeak), and hence in decreasing order of TFR nor¬ 
malization for fixed values of the other parameters. Thus 
models higher in the table require stronger selection and 
lower jy to achieve consistency with the data. 

• There are moderate degeneracies between the selec¬ 
tion factor and AM scatter (which influence the TFR 
scatter) and between u, j and the selection factor (which 
influence the normalization). The slope is the least con¬ 
straining of the three TFR characteriztics; it is satisfacto¬ 
rily reproduced in a large proportion of parameter space 
because the stellar mass range probed by the POT sample 
is relatively small (2.2 decades). 

• The maximum-likelihood values generated by each 
of the models are very similar, and correspond to a very 
good fit to the POT data. This demonstrates that our 
framework has sufficient flexibility to match the observed 
stellar mass TFR regardless of the AM scheme employed. 


5.2 Mass—size relation 

On constructing mock MSTs in our framework, we imme¬ 
diately find that the model is not capable of reproducing 
the observed intrinsic scatter, O.ITiO.OI dex. In particu¬ 
lar, the predicted scatter is always significantly larger (at 
least 0.28 dex), even when AM scatter is switched off and 
the selection factor maximized. The reason for this dis¬ 
crepancy is easily identified: in our model, in which the 
specific angular momentum of a galaxy is proportional to 
that of the halo in which it lives, disc size is proportional 
to halo spin. Thus the scatter in radius at fixed stellar 
mass receives a contribution equal to the scatter in halo 
spin at fixed AM parameter, which is roughly 0.27 dex 



0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.3 

(a) Slope 



(b) Intercept 



0 0.02 0.04 0.06 0.08 0.1 0.12 


(c) Intrinsic Scatter 


Figure 2. The distribution of slope, intercept and scatter 
values of 1000 mock TFR data sets, generated using the 
maximum-likelihood parameters of the Fpeak model (Table [sl, 
compared to the values for the observational data (red linS). 
There is good agreement for each power-law parameter. 


and therefore by itself exceeds the measured 
postpone further discussion of this issue to 


We use the ‘Peebles’ spi n (|Peebles|l969t for compatibility 
with previous work such as|Mo et al.| ( |I998| >. The alternative 
‘Bullock’ spin ( [Bullock et al.||200I| > has a typical scatter of 
around 0.29 dex and so would slightly exacerbate the discrep¬ 
ancy if used instead. 


value pq We 
Section (g^) 
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AM parameter 

AM scatter / dex 

3 

u 

Selection factor 

f^Deak 

< 0.18 (0.36) 

> 0.87 (0.63) 

n 4+0 -i(+0.3) 

> 1.09 (0.476) 

Pace 

< 0.16 (0.34) 

> 0.97 (0.67) 

+0.2(+0.4) 

- 

Pmax 

< 0.18 (0.36) 

> 1.0 (0.73) 

d d 

d d 

+ 1 

c5 

1 

- 

Mpeak 

< 0.18 (0.36) 

> 1.0 (0.72) 

+0.3(+0.5) 

^•^-0.2(-0.5) 

- 

M^acc 

< 0.20 (0.36) 

> 1.0 (0.77) 

+0.3(+0.6) 

^•^_0.2(-0.4) 

- 

Mnow 

< 0.24 (0.40) 

> 1.1 (0.81) 

fi 1 +0.2(4-0.6) 
^•-L_o.3(-0.5) 

- 


Table 2. Table of parameter constraints derived by comparing the model to the POT TFR. Numbers outside of brackets are Icr 
limits and numbers within brackets are 2(7 limits. Lower rows correspond to lower average halo concentration at given stellar mass. 
The constraints on both j and the selection factor depend sensitively on the prior range allowed for j (here 0 < j < 1.4); when 
j is allowed to take larger values the selection factor may be smaller. We quote limits on the selection factor only in cases where 
random selection (selection factor = 0.464) is ruled out at the > 2(7 level. 


AM parameter 

AM scatter / dex 

3 

u 

Selection factor 

P value 

Ppeak 

0.0 

0.86 

-0.4 

1.34 

- 0.9 

-Mnow 

0.4 

0.92 

0.2 

1.51 

- 0.9 


Table 3. Maximum-likelihood parameter values for the TFR and associated goodness of fit, for the Rpeak and Mnow models. ‘P 
value’ is roughly the fraction of mock data sets with {slope, intercept, scatter} values at least as far from the centroid as those of 
the data. All AM matching proxies are capable of generating good fits to the observed TFR. 



V 


Figure 3. Comparison of the u posterior probability distribu¬ 
tions (determining the halo response to disc formation) derived 
from the TFR for the Rpeak and Mnow models. Models that as¬ 
sign galaxies to haloes based on mass are consistent with mild 
halo contraction {u > 0) or expansion {u < 0); models that 
assign galaxies to haloes based on velocity require expansion. 


and in the remainder of this section investigate the ex¬ 
tent to which our model parameters can be constrained 
using only the slope and normalization of the MSR. 


To this end, we repeat the MCMC analysis described 
for the TFR in Section but replacing a galaxy’s Vso 
value by its scalelength, Pd, and calculating the likeli¬ 
hood for a particular set of parameter values by compar¬ 
ing the slopes and intercepts of power-law fits to mock 
data generated using the model to those of the data. We 
find that only j, the ratio of the specific angular momen¬ 
tum of the disc and halo, can be significantly constrained 
using this information, since this is mainly responsible 
for setting median disc size at fixed stellar mass, f values 
around 0.8 are favourecf^ (see Table and Fig. Q, sug¬ 
gesting that a small fraction of discs’ angular momentum 
may be transferred to their haloes. (Note that these con¬ 
straints on j include full marginalization over the other 
parameters in the model.) Slightly higher j values are 
preferred for velocity-based AM models than mass-based 
ones since the former tend to put relevant galaxies in 
more concentrated haloes, making those galaxies smaller 
at hxed angular momentum. In Appendix[A| (Fig. |A2| , we 
show and explain the effect of each of the model parame¬ 
ters on the predicted MSR. Fig. [^demonstrates that the 
model can ht well the slope and intercept of the relation 

Differences between the Peebles and Bullock spin param¬ 
eters are no more than ~ 10 per cent over the mass range 
of interest, so the specific definition of halo spin affects our 
constraints on j at the <10 per cent level. 


© 0000 RAS, MNRAS 000 , 000-000 
























12 H. Desmond and R. Wechsler 


Ypeak ' _, 

l^now ' ' 



1.2 1.4 


Figure 4. Comparison of the posterior probability distribu¬ 
tions of j (ratio of specific angular momentum of galaxy and 
halo) derived from the MSR, for the Rpeak and Mnow mod¬ 
els. Velocity-based AM models allow somewhat higher j than 
mass-based models. 


(using the best-fitting parameter values listed in the first 
row of Table , and illustrates its failure at reproducing 
the intrinsic scatter. 

Our choice to ignore the MSR scatter when deriving 
constraints on the model parameters corresponds to the 
assumption that a mechanism employed to resolve the 
discrepancy (see Section 6.2) would not significantly af¬ 
fect the median relation. We caution, however, that the 
constraints on j described above may be subject to sys¬ 
tematic uncertainties associated with our use of a model 
that does not adequately capture all aspects of the MSR. 
Furthermore, the constraint on j depends on the priors 
placed on the other model parameters. For example, re¬ 
stricting D and the AM scatter to small values and the 
selection factor to large values - in order to minimize the 
MSR scatter as much as the framework allows - tends to 
lower the preferred value of j, and reduce its statistical 
uncertainty. 


5.3 Correlation of Velocity and Radius 
Residuals 

If a galaxy’s stellar mass contributes significantly to its 
rotation velocity, one would expect smaller galaxies at 
fixed stellar mass to have larger Vso values, and hence 
that galaxies with different surface brightnesses would 
lie on systematically offset TFRs. In other words, there 
would be a significant anticorrelation between velocity 
and radius residuals at fixed stellar mass. This expecta¬ 


tion is not confirmed by the data (IZwaan et al. 1995 


McGaugh de Blok|1998||McGaugh|2005a||Courteau et| 

al.|2007| |. Here we explore this relation in the context of 
our model to elucidate its theoretical significance. 

We begin by taking the maximum-likelihood param- 


AM parameter 

3 

Fpeak 

n 0/1 +0-19(+0.38) 

Face 

+0.16(+0.27) 

'^•°°-0.15(-0.28) 

Frnax 

^ ^„+0.12(+0.25) 
^•®'-0.17(-0.26) 

Mpeak 

„.+0.18(+0.25) 

' -0.07(-0.16) 

Afacc 

. „„+0.14(+0.22) 

' -0.10(-0.19) 

Mnow 

n 09(4-0.19) 

^•'°-0.12(-0.19) 


Table 4. Table of constraints on j derived by comparing the 
model to the slope and intercept of the POT MSR. Numbers 
outside of brackets are Icr limits and numbers within brackets 
are 2(7 limits. 


eter values from Section |5.1| for a given AM proxy and 
generating 1000 mock data sets comprising radii and ve¬ 
locities of randomly chosen model galaxies with the same 
masses and velocity and radius uncertainties as the POT 
sample. We then fit the TFR and MSR of each mock data 
set with separate power laws, and calculate the residual of 
each galaxy as the difference between its log(V8o/km s“^) 
or log(Rd/kpc) value and the expected value for this 
quantity given its mass and the power-law fit. We then 
plot the residuals against each other and display the re¬ 
sulting contour plot (overlaying the points from ah 1000 
mock data sets) in Figs |6(aj] and [7(a)j for the Ipeak and 
Mnow models respectively. We quantify the strength of 
the correlation by calculating the Spearman’s rank cor¬ 
relation coefficient p for each mock data set and for the 
POT data, which specifies how close the two variables are 
to being monotonicahy related to each other, p — 1 cor¬ 
responds to a perfectly monotonic correlation, p = 0 to 
no correlation, and p = — 1 to a perfectly monotonic an- 
ticorrelation. 

Figs |6(c)| and |T(c)| (green histograms) compare the 
distributions of the 1000 coefficients thereby obtained to 
the corresponding value for the real data. Table (first 
column) displays the significance levels at which the hy¬ 
pothesis that the POT data was drawn from the model 
populations can be rejected, according to this test, for 
each AM proxy. We see that the parameter values that 
give the best fit to the TFR generate at least 8(7 ten¬ 
sion with the data, showing that our framework fails 
to capture observed velocity-radius correlations!^ The 
velocity-based AM models tend to produce a stronger 


The fact that the predicted radius residuals are signifi¬ 
cantly larger than those observed is a manifestion of the MSR 
scatter problem described in Section [5.2| Since p is calculated 
only from the ranks of the variables and not their absolute 
values, it would be unchanged if all radius discrepancies were 
scaled down by the same factor. 
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AM scatter / dex 

3 

u 

Selection factor 

P value (w/o scatter) 

P value (with scatter) 

l^Deak 

0.34 

0.60 

-0.2 

1.55 

> 0.9 

~ 10-13 

Afnow 

0.24 

0.5T 

-0.4 

1.55 

> 0.9 

~ 10-1° 


Table 5. Maximum-likelihood parameter values for the MSR aud associated gooduess of ht, for the Rpeak aud Mnow models. ‘P 
Value’ is roughly the fractiou of mock data sets with {slope, iutercept} (‘w/o scatter’) or {slope, iutercept, scatter} (‘with scatter’) 
values at least as extreme as those of the data. Our models are capable of reproduciug the slope aud iutercept of the observed 
MSR, but uot its small scatter. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

(a) Slope 


T5 


’i 


pfi 


B 

9 

z 



0.35 0.4 0.45 0.5 0.55 0.6 


(b) Iutercept 



0.1 0.15 0.2 0.25 0.3 0.35 0.4 

(c) lutriusic Scatter 


Figure 5. Same as Fig. but for the MSR iustead of the 
TFR (see Table |^. The slope aud iutercept agree well, but 
the predicted scatter is siguihcautly larger thau that of the 
data due to the relatively large scatter iu halo spiu. 



TFR params 

Optimum params 

hpeak 

12 

4.T 

Face 

14 

4.8 

Fmax 

IT 

5.1 

Mpeak 

8.1 

5.6 

Afacc 

8.2 

5.6 

ATnow 

8.4 

5.6 


Table 6. Siguihcauce level (iu a) at which the value of the cor- 
relatiou coefficieut p betweeu the radius aud velocity residuals 
of the POT data is iucousisteut with the theoretical predictiou. 
The hrst columu uses the parameter values that provide the 
best ht to the POT TFR; the secoud uses the parameter values 
that yield the weakest auticorrelatiou betweeu the radius aud 
velocity residuals aud heuce miuimize the siguihcauce of the 
disagreemeut with the data. 


anticorrelation and are therefore disfavoured by this test. 
The predicted value of p is largely driven by the tail of 
galaxies assigned to haloes with relatively high concen¬ 
tration (hence large Vso) and low spin (hence low Pd)- 
This tail contains ~ 10 per cent of the model galaxies. 

We note parenthetically that the discrepancy is even 
more pronounced for galaxies near the peak of the stel¬ 
lar mass-halo mass relation (M^ lO^^M©), because 

at these masses the galaxy’s relative contribution to the 
rotation velocity is maximized and the predicted anticor¬ 
relation of the residuals therefore strongest. Due to the 
paucity of POT galaxies near this mass, however, we are 
not able to strengthen our statistical conclusions by fo¬ 
cusing on this regime. 

We now investigate the possibility of resolving this 
discrepancy within our framework by adopting different 
parameter values. To this end, we repeat the MCMC 
analysis of Section]^ but replace the likelihood function 
with one determined by the residual correlation as op¬ 
posed to the TFR or MSR individually. In particular, we 
fit a normalized ID Gaussian to the distribution of p val¬ 
ues of our mock data sets, and take the likelihood to be 
the value of this Gaussian at the value of the real data. 

1 ( ax - 0 . 0485 )^ 

P(jp\d) oc —e (10) 

<7 

where p is the vector of model parameters, d is the p value 
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(a) l^peak; TFR parameter values 


A log Vgo I 

(a) Mnow; TFR parameter values 




(b) RpeakI optimum parameter values 


A log Vgo I 

(b) Mnow; optimum parameter values 



(c) 

Figure 6. Comparison of the Spearman’s rank correlation 
coefficient (p) between the radius and velocity residuals in the 
data and in 1000 mock data sets drawn from the Rpeak model. 
In Figs |6(a')] and [^(b)[ the black points are the POT galaxies 
and the contours enclose 20, 40, 60, 80, 90, and 95 per cent of 
the model galaxies. Fig. |6(a)| uses the parameter values that 
fit the TFR best; Fig. |6(b)| uses the optimum parameter values 
that maximize the p values of the mock data. Fig. |6(c)] shows 
the distribution of p across the mock data sets using the TFR 
(green) and optimum (blue) parameter values, and compares 
them to the p value of the real data (0.0485; red line). The data 
exhibit a significantly weaker correlation than is predicted by 
the model. 



(c) 

Figure 7. Same as Fig. but for the Mnow model. 
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of the POT data (0.0485), a is the standard deviation of 
a Gaussian fitted to the p values of the mock data sets, 
and p its mean. As in Section [5.11 we cap j at 1.4. 

For any AM parameter, we find that the predicted 
anticorrelation is minimized by maximizing j and the AM 
scatter, and minimizing ly. The discrepancy can be re¬ 
duced to the ~ 5cr level in this way [see Ta ble (sec¬ 
ond column), Figs |6(b)[ |7(b)| and Figs 6(c)[ 7(c) | (blue 
histograms)] at the price of predicting too low a TFR 
intercept and slope, and too high a scatter (P value of 
TFR fit < 10“^). In particular, almost any mock data set 
drawn randomly from our theoretical population would 
contain at least one galaxy with large positive velocity 
residual and large negative radius residual (correspond¬ 
ing to low halo spin and high concentration relative to 
the mean at that stellar mass); that the real data do not 
do so is therefore highly surprising. Further discussion of 
this problem, including its relation to previous work and 
potential solutions, may be found in Section [6.3| 


5.4 Dependence on Environment 

Haloes in different large-scale environments have differ¬ 
ent mass accretion histories (e .g. [Gottlober, Klypin “fc 


Kravtso^|200H |Wechsler et al.||2006| [Maulbetsch et al. 


2007), and galaxies with close neighbours may experience 
accelerated star formation due to tidal interactions (e.g. 


Kennicut^et^^y[1987jJMihos^&Ilernquis^|^^^^^^^^^^ 

&: Byrd||1996| I Alonso et al.||2004| [Barton et al.||2007| 


Behroozi et al. |2015| ). In addition, the relationship be¬ 
tween halo properties and galaxy star formation rates 
evolves differently for centrals and satellites. One might 
therefore expect galaxies in different environments to lie 
on systematically different TFRs. However, several previ¬ 
ous studies have found no evidence for such an effect. For 


example. Masters et al. (2006) find no difference between 


the TFRs of galaxies lying in 31 groups and clusters of 
varying sizes, and identify no correlation between TFR 
residuals and projected distance from the cluster centre. 


The analysis of Mocz et al. (2012) complements this work 


by splitting a large sample of SDSS galaxies into four 
bins of projected surface density and investigating their 
TFRs separately, finding no statistically significant dif¬ 
ference attributable to environmental effects. The data 


set of Gourteau et al. (2007) contains larger samples of 


both cluster and field spirals, and does not exhibit a sig¬ 
nificant difference between their respective TFRsHere 
we do not explicitly test for environment dependence in 
the POT data but instead ask whether our model pre¬ 
dicts a measurable dependence, and hence determine the 
information this relation could provide. 


However, there is some evidence for a systematic offset in 
the rotation velocities of satellite dwarf galaxies and isolated 
dwarf galaxies at low stellar mass jBlanton, Geha & West| 
200 ^. 


Rather than investigating environmental trends 
generically, we use as a measure of a galaxy’s environment 
whether or not it lies within the virial radius of another 
halo (i.e. whether it is a ‘satellite’ galaxy that resides in 
a ‘subhalo’, or a ‘central’ galaxy that resides in a distinct 
halo). For a given matching parameter, the standard AM 
approach makes a specific assumption about the galaxy 
growth for centrals versus satellites; we test the extent to 
which this is measurably manifest in the TFR. We begin 
by separating our model galaxies into satellites and cen¬ 
trals by splitting the halo population into subhaloes and 
main haloes as determined by ROCKSTAR. We adopt the 
maximum-likelihood parameter values from Section |5.1| 
and generate 1000 mock data sets, each with 156 galax¬ 
ies at the stellar masses of the POT galaxies. We then fit 
a separate power law to the TFRs of the satellite and 
central galaxies of a given mock data set, and record the 
distributions of the differences of their slope, intercept 
and scatter values. Finally, we calculate the significance 
levels at which no difference between the populations can 
be excluded by the model; in other words, we calculate 
the confidence with which one would reject the hypoth¬ 
esis that the real data was drawn from the theoretical 
populations if one were to perform this analysis on the 
POT data and find no difference between the satellite and 
central TFRs. 

Fig. shows the separate theoretical TFRs for the 
best-fitting Hpeak and Mnow models, and Table quan¬ 
tifies their differences (along with the other AM mod¬ 
els). We find that although satellite galaxies have lower 
halo Mnax than centrals of the same stellar mass (due to 
stripping of their outer regions during infall), they tend 
to have slightly higher Vso values. This is because the 
rotation curves of subhaloes peak at smaller radii than 
distinct haloes (due to their higher concentration), closer 
to the Rso values of the galaxies they host. In addition, 
satellite galaxies are predicted to exhibit a larger scat¬ 
ter in Vso, since the strength of stripping varies between 
subhaloes. This effect is larger for models that match to 
halo velocity than mass because these cause more satel¬ 
lite galaxies’ stellar masses to be set before stripping. 
However, these differences are not significant: were one 
to perform this test on the real data and find no differ¬ 
ence between the two TFRs (in particular their intrinsic 
scatters), the most one could say is that the 14cc and 
Hpeak models are disfavoured at the 2cr level. 

The main reason for the statistical insignificance of 
the difference between the theoretical satellite and cen¬ 
tral TFRs is the size of the sample: with satellite frac¬ 
tions around 10-35 per cent in this stellar mass range, 
the satellite-only mock data sets are small, resulting in 
a large variation among the parameters of the power- 
law fit between them. However, the framework we have 
set up is easily generalized to permit the analysis of a 
hypothetical data set with similar characteriztics to the 
POT sample but more galaxies, allowing us to estimate 
the number of galaxies needed to definitively manifest a 
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■oglMstar^ sun! 


(a) V^eak 



(b) Mnow 

Figure 8. Comparison of the predicted satellite-only and 
central-only TFRs, using the best-fitting parameter values as¬ 
sociated with the Rpeak and Mnow models respectively. The 
points show the median velocities of the model galaxies in 
1000 mock data sets, and the error bars show the interquartile 
range. TFRs for satellite galaxies only are predicted to have 
somewhat higher normalizations and intrinsic scatters than 
those for central galaxies only. 


nificant difference remains in the scatter. The separation 
expected for the Rpeak model is over 5cr. Were such a 
homogeneous data set available which used assumptions 
and definitions consistent with our model, an analysis of 
the present type would yield useful information about 
the general validity of our framework, and might pro¬ 
vide another way to distinguish between the various AM 
schemes. 

We caution that our methodology assumes that the 
AM process itself makes no distinction between the satel¬ 
lite and central populations, so that they can be cleanly 
separated after the matching has been performed. It is 
also possible that satellite galaxies should instead be as¬ 
signed to different haloes at the outset, or that effects not 
taken into account by AM (e.g. significant continued star 
formation after accretion, accelerated star formation due 
to galaxy-galaxy interactions, or systematic variations in 
halo shape with environment) are significant. In this case 
our model would need to be generalized for agreement 
with the data to be expected. 


6 DISCUSSION 

We have constructed a general framework for assessing 
whether AM-based galaxy formation models are capa¬ 
ble of reproducing the detailed characteriztics of the ob¬ 
served Tully-Fisher and mass-size relations, and deter¬ 
mined what information can be gleaned from them. In 
Section |6.1| we compare our parameter constraints to 
those obtained by studies in the literature. We then pro¬ 
vide in Sections |6.2| and |6.3| respectively, a detailed dis¬ 
cussion of the two major problems that our analysis has 
revealed: the theoretical MSR scatter is too large, and the 
velocity and radius residuals are too strongly anticorre¬ 
lated. In particular, we show that the solution to both 
of these problems (in addition to several others already 
known) requires a correlation between stellar surface den¬ 
sity and enclosed dark matter mass that is not present in 
models currently under investigation. 

6.1 Parameter Constraints and Comparison 
with the Literature 


difference between the satellite and central TFRs. In Ta- 


The simplest AM models (which assume adiabatic or 
near-adiabatic contraction and do not take account of 
selection effects) give galaxies velocities that are too high 
for their stellar mass. This result is present in several 


set comprising 700 galaxies with log(M*/MQ) uniformly 

other TFR studies (e.g. Dutton et al. 2007 

Trujillo- 

distributed in the range of the POT sample (9-11.2) and 

Gomez et al. 201 1| |Desmond 2012 McGaugh 

2012b), 


with log(V80/km s“ ) uncertainties equal to the POT RMS 
value (0.0346). In this case, we would expect to measure 
a significant (> 3cr) difference between the two TFRs 
for the velocity-based AM models, but not for the mass- 
based ones. We note that although there is a difference 
in the mean relations between the two populations, it is 
still too small to be robustly distinguished; the most sig- 


and its quantitative solution likely requires not only that 
haloes contract less strongly than prescribed by the adi¬ 
abatic formalism (in agreement with Dutton et al.|2007 ), 
but also that observational TFR studies select galaxies 
with systematically low velocities for their stellar mass. In 
other words, rotation-supported spirals must live in low- 
mass or low-concentration haloes at fixed stellar mass. 
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This trend has been suggested by other theoretical con¬ 
siderations ( Rodriguez-Puebla et al.|2011 ) in addition to 
measurements of satellite kinematics (Wojtak & Mamon 


2013), but our model allows us to investigate quantita¬ 


tively the required magnitude of the effect and its de¬ 
pendence on other galaxy formation parameters. In the 
context of the TFR, the required magnitude of selection 
is degenerate with the ratio of the specihc angular mo¬ 
mentum of the disc and halo. Improving our knowledge of 
j would therefore afford a tighter constraint on the selec¬ 
tion factor. We have checked that the difference between 
the haloes of ‘spirals’ and ‘other types’ suggested by our 
analysis is not in obvious disagreement with the measure¬ 
ments of Wojtak & Mamon (2013), but it remains to be 


seen whether the required magnitude of the difference be¬ 
tween the selected subsample and the full population is 
theoretically reasonable and can be corroborated in detail 
by independent measurements. 

Reproducing the small scatter in the stellar mass 
TFR requires that little scatter is associated with AM 
itself - that is, the stellar mass is in near-monotonic cor¬ 
respondence with halo mass or velocity. In particular, the 
scatter in this relation must be less than around 0.35 dex 
(2cr) for all models tested. This is consistent with the re¬ 
sults of Reddick et al. (2013), who find that an AM scat¬ 
ter of ~ 0.2 dex is required for agreement with satellite 
fraction and galaxy clustering measurements. We have 
also tried focusing attention on the least massive POT 
galaxies in order to set a limit on the AM scatter that 
applies specihcally at low stellar mass, where constraints 
from alternative methods are weak. We found that the 2(7 
upper limit on the scatter is around 0.7 dex for a subset 
of the POT data with median log(M*/M0) of 9.6. How¬ 
ever, this constraint is weaker than that obtained from 
the full data set only because fewer galaxies are used, 
and we are unable to probe even lower stellar masses 
due to the paucity of galaxies in the POT sample with 
log(M*/Mo) < 9.5. Applied to a future TFR data set 
with more galaxies, however, this method should clearly 
reveal whether AM scatter values required for clustering 
and the TFR are consistent, and set robust constraints 
on the scatter at low stellar mass. Constraining the scat¬ 
ter at low mass is difficult with clustering because the 
bias function is not a strong function of halo mass in 
this regime, and from group statistics due to the scarcity 
of samples of robustly identified groups with centrals of 
Mvir < 10^^ Mq. 

Turning now to the MSR, we see from Table that 
the ratio of the specific angular momentum of a disc and 
its halo [j) must be in the approximate range 0.6-1.2 
to ensure the correct normalization of the relation. This 
result suggests that the baryonic material retains most 
of its primordial angular momentum when it collapses 
into a disc, and further reinforces the need to resolve 
the overcooling problem that has plagued hydrodynami- 
cal galaxy simulations ( Navarro fc Whit^|I994 Navarro 


Steinmetz||I997| |Dekel Maller||2Q03| |Scannapieco et 
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al.|20I2 ). Our result is in rough agreement with Kravtsov 


(2013) who use a simplihed version of the Mo et al. (1998) 
method to compare the theoretical and observed sizes 
of galaxies with a range of morphologies across a wide 
range of stellar mass. We anticipate that our constraint 
will provide a useful touchstone for prescriptions of bary¬ 
onic dissipation and infall implemented in hydrodynam- 
ical galaxy simulations. 


6.2 The Excess Intrinsic Scatter Predicted for 
the Mass—Size Relation 


A signihcant problem identihed in this work (shown most 
clearly in Fig. 5(c)| is that our model predicts that the 
MSR scatter should be signihcantly larger than is ob¬ 
served. The smallest scatter we are able to obtain is 
around 0.28 dex (by matching to Mnow and setting AM 
scatter = 0, selection factor = 7r/2), while that of the 
POT data is 0.17 ± O.OI dex. The theoretical scatter in 
disc scalelength Rd at fixed stellar mass is set by five 
factors, which are assumed to be uncorrelated. 


(i) The scatter in halo spin A at hxed virial mass (at 
fixed j, A sets a galaxy’s angular momentum and hence 
regulates its size). 

(ii) The scatter in halo concentration at hxed virial 
mass (more concentrated haloes generate larger gravita¬ 
tional potential gradients in the inner regions where the 
disc is situated, boosting the galaxy’s rotation velocity 
and hence reducing its radius at hxed angular momen¬ 
tum). 

(hi) The scatter in present-day halo mass at hxed stel¬ 
lar mass (set by the AM prescription). 

(iv) The scatter in bulge mass fractions. 

(v) The impact of selection effects. 


Of these, the third can be eliminated by adopting the 
Mnow am model with zero scatter, the fourth is small, 
and the hrst contributes a hxed 0.27 dex (in our model, 
Rd (X X). This is already more than is allowed for by the 
data. When galaxies are selected on the basis of their Vso 
values, as we have assumed, selection effects are capa¬ 
ble of removing around 0.21 dex in quadrature. Scatter 
in concentration largely makes up the remaining ~ 0.22 
dex. These results agree with the broad-brush analysis 
of Kravtsov (2013), who hnds the observed relation to 
have a scatter of 0.2 dex (the difference to POT’s 0.17 dex 
being largely attributable to the subtraction of observa¬ 
tional uncertainties in the latter), but does not include 
the additional theoretical scatter expected from concen¬ 
tration, stellar mass fraction or bulge mass fraction, and 
does not consider selection effect srM 


[Courteau et al.| ( |2007| ) hnd an even smaller MSR scatter 
of ~ 0.14 dex, although the use of a different htting method, 
and luminosity as opposed to stellar mass, make this result 
not directly comparable to ours. 











































18 


H. Desmond and R. Wechsler 



Slope 

Intercept 

Scatter 

Vpeak 

0.23 

1.6 

2.3 

Vacc 

0.29 

1.5 

2.0 

Vinax 

0.12 

1.2 

1.3 

Mpeak 

0.48 

1.1 

0.95 

Mace 

0.34 

1.0 

0.44 

Afnow 

0.10 

0.82 

0.61 


Table 7. Table of significance levels (in a) at which the 
satellite-only and central-only TFRs of the POT data set are 
expected to be distinguishable, according to each power- 
law parameter, using the parameter values that provide 
the best fit to the POT TFR. AM models that match to 
halo velocity predict a larger difference between the two 
TFRs than those that match to halo mass, although the 
difference is not significant in any case. 



Slope 

Intercept 

Scatter 

f^eak 

1.9 

1.8 

5.1 

Vacc 

1.8 

l.T 

4.2 

Vrnax 

1.0 

1.3 

3.1 

Mpeak 

1.4 

1.0 

2.2 

Afacc 

1.3 

1.1 

1.4 

AJnow 

0.025 

1.0 

0.41 


Table 8. Same as Table but for a hypothetical data 
set comprising TOO galaxies with stellar masses uniformly 
distributed in the range 9 < log(M*/MQ) < 11.2. In this 
case, the satellite and central-only TFRs are expected to 
be distinguishable on the basis of their intrinsic scatters (at 
the > 3(7 level) for all of the velocity-based AM models. 


This result implies that the simplest angular 
momentum-based methods for assigning disc sizes do not 
agree with the observed scatter in the MSR. This prob¬ 


lem was noted by de Jong Sz Lacey (2000) and further 


quantified by Gnedin et al. (200T), who noted that the 


width of the A distribution of discs can be no larger than 
~ O.IT dex. The discrepancy would be exacerbated by 
introducing scatter into the ratio of the specific angular 
momentum of the disc and halo (presumably expected at 
some level due to differences in their formation). We now 
describe three specific ways in which this problem could 
be resolved whilst retaining the notion that galaxy sizes 
are set by angular momentum partition between the disc 
and halo. 

The first is to argue that selection effects in the data 
play a more significant role than is allowed for in our 
model. In particular, we could reduce the predicted MSR 
scatter by more than 0.21 dex if we assume that the quan¬ 
tities by which the POT galaxies were chosen for inclusion 
are correlated more strongly with galaxy size or host halo 
spin than with Vso- For example, it has been suggested 
that low angular momentum discs are more prone to the 
development of a bar instability, and may therefore be¬ 
come early-type galaxies (e.g de Jong 1995| Dalcanton 
Spergel fc Summers|199T Mo et al.|1998 ). Three points 
however, argue against this as a solution to the problem 
First, we have tried constructing a simple toy model in 
which our model galaxies are selected not on the basis of 
Vso alone but rather a linear combination of Vso and Rd, 
and have not been able to reduce the predicted MSR scat¬ 
ter to the observed value while preserving agreement with 
the TFR. It is possible to lower the MSR scatter without 
affecting that of the TFR by preferentially eliminating 
low-spin haloes at fixed and Vso; however, conditional 
selection of this type is somewhat fine-tuned and we do 
not investigate it further here. A more physically moti¬ 
vated prescription that gives similar correlations could be 
interesting to investigate in future work. Secondly, a res¬ 
olution of the MSR scatter problem in terms of selection 


effects would imply that a galaxy sample chosen using a 
broader set of selection criteria would exhibit significantly 
larger scatter than the POT sample. As noted above, how¬ 
ever, the Kravtsov (2013) sample includes galaxies with a 
range of morphologies yet does not find an intrinsic MSR 
scatter discrepant with that of POT. Finally, reducing the 
~ 0.35 dex predicted scatter to the required O.IT dex by 
selecting only the ~ 50 per cent with largest radii would 
impart a strong skew on the resulting distribution, which 
would therefore have a tail towards high radius. This is 
not apparent in the observational datap^l 

Second, one might attempt to resolve this problem 
by using the fact that star formation proceeds faster in 


higher-density regions (Schmidt 


1959 


Kennicutt 


1989). 


This suggests that at fixed baryonic mass, larger discs 
would have smaller stellar masses, which could reduce the 
width of the size distribution at fixed stellar mass rela¬ 
tive to the distribution at fixed total baryonic mass. This 
mechanism could therefore be tested by comparing the 
scatters of the and Mbaryon-Rd relations. How¬ 

ever, it cannot be accommodated within the basic AM 
framework where stellar mass is set purely by halo mass 
or velocity. 

A final way of lowering the predicted scatter in the 
MSR is to postulate a correlation between a galaxy’s stel¬ 
lar surface mass density (or equivalently its size at fixed 
stellar mass) and the mass or concentration of the halo in 
which it lives. Suppose that at fixed stellar mass, larger- 
than-average galaxies are made to live in more massive 
haloes. The effect of this would be to increase the ro¬ 
tational velocity of such galaxies. At fixed disc angular 
momentum, this would reduce those galaxies’ radii, mov¬ 
ing them closer to the average. Thus the magnitude of 


This argument applies also to the TFR, and suggests that 
future data sets with many more galaxies than POT may be 
able to constrain the impact of selection by means of the third 
moments of the velocity and radius distributions in bins of 
stellar mass. 
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radius discrepancies would be decreased. An equivalent 
way to view this is to suppose that at hxed halo mass, 
larger galaxies are made to have larger stellar mass. This 
would push a larger-than-average galaxy to the right on 
the plane, reducing the difference between its size 

and the average for that stellar mass. From either per¬ 
spective, we see that the intrinsic scatter in the MSR 
would be reduced. 

This is not the only observation which suggests a 
correlation between surface mass density and enclosed 
dark matter mass, and evidence for such a correlation 
exists for early as well as late-type galaxies. It has been 


shown by Sonnenfeld et al. (2015) using weak leasing 


and kinematical observations of ellipticals that the in¬ 
ferred dark matter mass within the central 5 kpc anticor¬ 
relates with surface brightness. In addition, three very 
massive and high surface brightness ellipticals have been 
observed to have velocity dispersions that fall off in a Ke- 
plerian fashion from the centre of the galaxy, indicating 
that the enclosed dark matter is dynamically insignih- 
cant ( Romanowsky et al.|2003 ). The correlation has been 
documented for a wide range of galaxy types in |Gen-| 


tile et al. (2009), and we show in Section 6.3 that the 


observed relation between the residuals of the TFR and 
MSR is yet another manifestation of this effect. The most 
pristine form of the correlation, however, is the ‘mass 
discrepancy-acceleration relation’ for spirals, which sug¬ 
gests a tight correlation between baryonic surface mass 
density and enclosed dark matter mass (|Sanders| ~990| 


McGaugh|[2^ 

Gaughliobf"^^ 


[Tiret fc Combes||2009| [Famaey &: Mc-| 


We identify two mechanisms whereby a relation of 
this form could be introduced into our model. The first 
is to postulate a correlation between the AM scatter (the 
scatter in stellar masses of galaxies associated with haloes 
of fixed mass or velocity) and halo spin, which is the main 
determinant of galaxy size. In particular, we might sup¬ 
pose that the galaxies with the greatest stellar masses for 
their halo mass live in the most rapidly rotating haloes, 
tending to make them larger than average. In this case, 
increasing the AM scatter would decrease the scatter in 
the MSR for the reason described above. We have inves¬ 
tigated this possibility using a simple toy model, which 
suggests that an AM scatter of around 0.3 dex may be 
required to bring the MSR scatter down to the observed 
level Fu This is inconsistent with the conditional SMF 


measurements of Reddick et al. (2013). Furthermore, a 


perfect correlation between AM scatter and halo spin is 
unlikely, and a noisy correlation would require an even 


Indeed, it is evidence for such a correlation that gives Mod¬ 
ified Newtonian Dynamics (MOND) its phenomenological ef¬ 
ficacy. 

Alternatively, one could anticorrelate abundance matching 
scatter with halo concentration. This would require a larger 
AM scatter to be effective because disc size at fixed stellar 
mass is affected less strongly by concentration than spin. 


larger AM scatter to be effective. Nevertheless, we can¬ 
not rule out this scenario as potentially offering a solution 
to this problem. 

An alternative way to enforce a correlation between 
stellar surface mass density and dark matter mass would 


be to perform ‘conditional abundance matching’ (Hearin 


& Watson 2013 Hearin et al. 2014). In this case, one 


would no longer assign disc sizes via angular momentum 
considerations, but instead use directly the observed dis¬ 
tribution of disc sizes as a function of stellar mass. The 
radii of model galaxies would be assigned in monotonic 
or near-monotonic correspondence with the present-day 
mass or concentration of the haloes in which they live. 
The disadvantage of this method is that it loses the the¬ 
oretical motivation of the Mo et al. (1998) model (that 
the baryons and dark matter share angular momentum), 
but, if successful, could provide an important clue to un¬ 
derstanding the galaxy-halo connection. 


6.3 The Excessive Anticorrelation between 
Radius and Velocity Residuals 

The second issue illuminated by our work is that the AM 
plus disc size models we consider predict too strong an an- 
ticorrelation between a galaxy’s radius discrepancy (the 
difference between its disc scalelength and that expected 
given its stellar mass and a power-law ht to the full MSR) 
and its velocity discrepancy (the analogue for the TFR). 
This is illustrated by Figs and We have demonstrated 
that models which match the TFR predict an anticorrela¬ 
tion too strong at the > 8cr level, and the discrepancy can 
be reduced (though never below ~ 5cr within the present 
framework) only by adopting parameter values inconsis¬ 
tent with the TFR alone. This suggests that models of 
this type, at least with standard assumptions about the 
galaxy-halo connection, are incapable of capturing in de¬ 
tail the correlations observed between galaxies’ sizes and 
rotation speeds. Here we compare our conclusion to oth¬ 
ers in the literature, and provide suggestions for how the 
discrepancy might be resolved. 

One of the best-known studies of the correlation of 


velocity and radius residuals is Gourteau & Rix (1998), 


in which it is argued that the observed level of correlation 
requires that the disc contributes (60 ± 10) per cent of 
the total rotation velocity at 2.2 disc scalelengths. In this 
case, according to their model, the increase in rotation 
velocity due to the baryons as the disc’s size decreases at 
hxed stellar mass is fully compensated for by a reduction 
in the amount of dark matter mass within R2.2. In our 
model the discs contribute on average an even smaller 
fraction of the total rotation velocity (due to the fact that 
AM induces stellar mass fractions far below the cosmic 
baryon fraction), yet result in a residual correlation that 
is in stronger disagreement with the data. We identify 
three reasons for this difference. 


First, the Gourteau & Rix (1998) model uses the fact 


that halo spin is uncorrelated with concentration to argue 
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that galaxy size is too. However, halo spin sets only the 
angular momentum of the disc, and does not fully specify 
its size; in more concentrated haloes, discs rotate faster 
and are therefore smaller at fixed angular momentum. 
This effect intensihes the predicted anticorrelation of the 
radius and velocity residuals: not only is the galaxy’s own 
contribution to the velocity larger when the galaxy is 
smaller, the halo’s is likely to be too. Secondly, a signif¬ 
icant fraction of our model galaxies are larger than the 
fiducial 3 kpc assumed by Courteau Rix| ( |1998| ) , and we 
measure the rotation velocity at Rso (about 3Rd) instead 
of i^2.2. This means that we are probing a larger physi¬ 
cal radius, often around the peak of the halo’s rotation 
curve, and hence the halo does not provide significantly 
greater rotation velocity for larger discs. Finally, the de¬ 
gree of residual anticorrelation found in the sample to 


which Courteau & Rix (1998) compare their results is 


significantly larger than ours (the slope is —0.18±0.05 in 
their ‘MAT’ sample). The statistically insignificant cor¬ 
relation observed in the POT sample appears more repre¬ 
sentative of TFR studies in general (|Verheiien|2001| |Mc- 


Gaugh 2005a| Courteau et al. 2007| ^yes et al.||2Qll|). 

Furthermore, our analysis improves upon |Courteau fc Rix| 
(1998) by modelling all correlations between disc and halo 


properties (stellar mass, disc scalelength, halo mass, con¬ 
centration, and spin), including their scatter, in a fully 
self-consistent and integrated way that relies only on the 
AM ansatz and the assumption that halo spin regulates 


galaxy size. We conclude that the Courteau & Rix (1998) 


solution to the residual correlation problem is not prac- 
ticableF^ 

Several other methods by which one might attempt 
to resolve the discrepancy can be seen to lead to or ex¬ 
acerbate conflict with other observations. For example, 
making discs larger at fixed stellar mass would reduce 
the galaxy’s contribution to Vso but generate disagree¬ 
ment with the normalization of the MSR. Increasing the 
dark matter mass in the central regions would aggra¬ 
vate the ‘cusp/core’ and, potentially, the ‘too big to fail’ 
problems. Finally, reducing stellar mass fractions would 


intensify the missing baryon problem (McGaugh 2008 
McGaugh et al.|2010 ). Furthermore, implementing these 


would likely require a break with the AM ansatz. Al¬ 
ternatively, one might hypothesize that selection criteria 
in typical TFR surveys conspire to reduce the correlation 


Other studies in the literature attempt to eliminate the 
residual anticorrelation by postulating and fine-tuning a cor¬ 
relation between stellar mass and baryonic surface density in¬ 
duced by a star formation threshold (e.g. [Firmani Sz AvilaT] 
|Reese|[2000[ [Dutton et al.|[^07| >. Such models nevertheless 
predict a non-zero anticorrelation of the residuals, so it is un¬ 
clear whether they are in statistical agreement with the data. 
Furthermore, they predict that baryonic TFR residuals would 
be significantly anticorrelated, in conflict with some observa¬ 
tions ( |McGaugh||2005^ . In any case, this mechanism is not 
available within the basic AM framework. 


by eliminating small galaxies with high velocity and large 
galaxies with low velocity. This is possible but not well 
motivated, and is not guaranteed to be successful if the 
scatter in the residual correlation is insufficiently large. 

The problem could be alleviated by imposing a cor¬ 
relation between stellar surface density and halo mass or 
concentration, as described in Section |6.2| As discussed 
above, in models like ours which set galaxy size by an¬ 
gular momentum considerations, smaller galaxies tend 
to reside in more concentrated haloes. However, the ob¬ 
served lack of correlation of the residuals indicates that 
at fixed stellar mass, smaller galaxies must reside in less 
massive or concentrated haloes so that the increase in ro¬ 
tation velocity due to the galaxy is offset by a decrease 
in that due to the halo. This conclusion was also reached 
by Gnedin et al. (2007), who found agreement with the 
observed correlation of the residuals when central disc- 
to-halo mass fraction rud and surface stellar mass den¬ 
sity E* = M* R^^ were related by ma oc E^'®^ (assuming 
standard adiabatic contraction). It remains to be seen 
whether either of the strategies described in Section [6.2| 
for introducing the required correlation within the AM 
framework (or indeed any other) are capable of simul¬ 
taneously resolving both problems identified here, repro¬ 
ducing observations of ellipticals, and receiving a priori 
theoretical support. 


7 SUGGESTIONS FOR FURTHER WORK 

We identify several ways in which this study could be 
usefully extended or its results used to further improve 
our knowledge of galaxy formation. 

• We have used the TFR and MSR to constrain the 
ratio of disc to halo specific angular momentum, the re¬ 
sponse of the halo to disc formation, the scatter in the 
relation between stellar mass and halo mass or rota¬ 
tion velocity, and the impact of selection criteria. These 
constraints require theoretical motivation from detailed 
hydrodynamical simulations of galaxy formation, and 
should be checked for consistency against other observa¬ 
tions (e.g. the Faber-Jackson relation and Fundamental 
Plane of ellipticals) using a similar framework. 

• Enlarging the observational TFR data set could sig¬ 
nificantly improve constraints on the scatter between 
galaxy and halo properties. At present the AM scatter 
is better constrained by clustering than by the TFR, but 
a larger TFR sample would allow a stronger test of con¬ 
sistency between the two methods. Extending the TER 
sample to low mass would enable the mass dependence 
of the scatter to be investigated in a regime that is not 
accessible to clustering measurements. An analysis of the 
baryonic TER along the lines described here would, when 
feasible, likely provide even more constraining power. 

• As described in Section [5^ a sample of around 700 
galaxies split into satellites and centrals should reveal an 
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significant difference between the TFRs of the two popu¬ 
lations if galaxies populate haloes according to a velocity- 
based AM scheme. Quantifying this difference would pro¬ 
vide a direct handle on potential differences between the 
galaxy-halo connections of centrals and satellites. 

• We have argued that several problems associated 
with galaxy phenomenology require a correlation between 
galaxy stellar surface mass density and halo mass or con¬ 
centration that is not induced naturally by simple one- 
parameter AM models. This could be produced either by 
correlating AM scatter with halo spin or concentration, 
or by using conditional abundance matching to correlate 
galaxy size with halo mass or concentration at hxed stel¬ 
lar mass. An obvious next step is to implement such mod¬ 
els to see whether they are capable of resolving the prob¬ 
lems described in Section and, if so, what strength of 
correlation is required. A correlation between disc scale- 
length and halo mass would suggest that galaxies of dif¬ 
ferent size but the same stellar mass would cluster differ¬ 
ently, which could be tested in SDSS data. A correlation 
between halo concentration and galaxy size would also 
likely imply some degree of galactic conformity in the 
theoretical population, which could be compared to ob¬ 
servations ( Weinmann et al.|2006 Hearin, Watson & van 


den Bosch|2Q15 Paranjape et al.|2Q15 ). It would then be 

incumbent upon hydrodynamic simulations or detailed 
analytic calculations to explain how the correct correla¬ 
tion between stellar surface mass density and dark matter 
mass comes about. 


8 SUMMARY AND CONCLUSIONS 

We have explored a range of plausible assumptions about 
the galaxy-halo connection, implementing AM models 
that have previously been shown to match a number of 
observed galaxy statistics in the local Universe. We in¬ 
vestigated the extent to which these models are able to 
reproduce the stellar mass Tully-Fisher and mass-size re¬ 
lations for local galaxies, using a homogeneously selected 
data sample from the SDSS. Our main hndings are as 
follows. 

• A model that matches galaxy stellar mass to max¬ 
imum halo velocity (either today, at accretion, or at 
the point in its merger history when this quantity was 
largest) can fit the normalization of the stellar mass TFR 
in the range 9 < log(M*/Mo) <11.2 only if the halo ex¬ 
pands in response to disc formation. (A small amount of 
contraction is permitted if Um and take their lowest 
values allowed by the Planck uncertainties.) For a model 
that matches instead to halo mass, no contraction or ex¬ 
pansion is favoured. In either case, ‘standard’ adiabatic 
contraction is ruled out at the > 4cr level. 

• The small intrinsic scatter of the stellar mass TFR 
places an upper limit on the AM scatter (< 0.35 dex at 
2(7) and suggests that galaxies selected for observational 
TFR studies preferentially populate the low end of the 


rotational velocity function at hxed stellar mass. Once 
these constraints are satished, an AM model that matches 
to any of the properties considered is able to successfully 
reproduce the observed TFR. 

• The slope and normalization of the MSR can be 
matched simultaneously with the TFR, and constrain the 
ratio of the specihc angular momenta of the baryonic and 
dark matter mass to be in the approximate range 0.6—1.2. 
However, the intrinsic scatter predicted for the MSR is 
signihcantly larger than that in the data even when AM 
scatter is switched off (at least 0.28 dex, compared to 0.17 
dex). This appears to be a fundamental consequence of 
setting galaxy size proportional to halo spin, which itself 
has significant scatter. 

• The observational data show no statistically signif¬ 
icant correlation between the residuals of the MSR and 
TFR. However, the model predicts an anticorrelation be¬ 
cause smaller galaxies at hxed stellar mass generate a 
larger rotation velocity. The Spearman’s rank correlation 
coefficient is found to be at least 8cr discrepant with the 
data for parameter values which generate a good ht to 
the TFR. This discrepancy can be reduced to around 5(7 
only by adopting parameter values that do not generate 
agreement with the TFR alone. 

• Both of the above problems (in addition to several 
others already known) imply a correlation between stel¬ 
lar surface mass density and enclosed dark matter mass 
which is stronger than that induced by the current stan¬ 
dard range of AM models. This could be introduced by 
correlating AM scatter ad hoc with halo spin or concen¬ 
tration, or by using ‘conditional abundance matching’ to 
connect galaxy size to halo mass or concentration at hxed 
stellar mass. It remains to be seen whether a model along 
either of these lines can be theoretically motivated and 
used to successfully resolve the discrepancies described 
here. 

• In the context of our model, splitting the data sam¬ 
ple we have investigated into satellite and central galax¬ 
ies would not be expected to yield distinguishable TFRs. 
However, we predict that a sample with uniform selec¬ 
tion and at least 700 galaxies could be resolved into 
separate satellite and central TFRs at the > 3cr level 
if galaxies populate haloes according to a velocity-based 
AM scheme. Analysis of a future data set in the way we 
describe would provide another crucial test of the AM 
framework, and allow us to further constrain its degrees 
of freedom. 
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APPENDIX A: MODEL FITS TO THE TER 
AND MSR 


We describe and explain here the effect of each of the 
individual model parameters on the TER (see Eig. A1) 
and MSR (see Eig. A2). 


• Reducing j increases the normalization and scatter 
of the TER, and decreases the normalization of the MSR. 
Lower j corresponds to smaller disc angular momentum, 
and hence smaller size for fixed values of other disc and 
halo parameters. This increases the disc’s contribution to 
the rotation curve, which contributes additional scatter 
to Vso according to the scatter in halo mass, concentra¬ 
tion, and spin, and stellar bulge-to-disc ratio, at fixed 
stellar mass. 

• Increasing v causes more dark matter mass to be 
pulled within Rso, increasing the normalization of the 
TER. This lowers slightly the normalization of the MSR, 
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since faster-rotating discs are smaller when other proper¬ 
ties are hxed, and increases its scatter. 

• Increasing the AM scatter generates additional vari¬ 
ation in halo mass and concentration at hxed stellar mass, 
increasing the scatter in the TFR and MSR. It also re¬ 
duces the small amount of curvature in the predicted 
TFR, and decreases the MSR slope. 

• Reducing the selection factor weakens the preference 
for low-Vso galaxies to be selected at hxed stellar mass. 
This increases the normalization of the TFR, reduces the 
normalization of the MSR (high-Vso galaxies tend to be 
small), and increases the scatter in both. 


APPENDIX B: ASSUMPTIONS AND 
SYSTEMATIC UNCERTAINTIES 

We review here the most signihcant assumptions in our 
analysis. We discuss the consequences of reasonable de¬ 
viations, and, where possible, provide quantitative esti¬ 
mates of the systematic errors that may be induced, in or¬ 
der to demonstrate the robustness of our results to them. 
Overall, we believe the systematic error to be under con¬ 
trol at the 10 per cent level or better. 

Fraction of galaxies removed by selection cuts', we have 
assumed that the exclusion from the POT sample of galax¬ 
ies with rotation curves poorly ht by an arctan function 
does not signihcantly bias the halo population in which 
these galaxies live, i.e. that the deviation of a galaxy’s 
rotation curve from an arctan function is not strongly 
correlated with halo properties. In reality, some galaxies 
likely failed this cut because they were early-type and 
did not show strong signs of rotation - as discussed in 
Section early-type galaxies are expected to reside in 
systematically different haloes than late-type galaxies at 
fixed galaxy mass. However, we do not expect this to af¬ 
fect our main conclusions since the fraction of galaxies 
removed as a function of mass is degenerate with our se¬ 
lection factor parameter. Reducing the fraction of galax¬ 
ies retained during the selection step would merely lower 
the preferred value of the selection factor. 

Modelling of selection effects', we have assumed that se¬ 
lection effects can be modelled reliably via a correlation 
between observed disc axis ratio h/a (the primary quan¬ 
tity used for selection in POT) and Vso- This is not incon¬ 
sistent with the more physical assumption of a correlation 
between a direct morphological proxy (e.g. Sersic index 
or true axis ratio) and any alternative measure of dynam¬ 
ical mass (e.g. Umax or halo Mvir), which would imply a 
correlation (albeit weakened) between h/a and Vso- Our 
methodology will therefore capture the first order effect 
of any such correlation. Although the impact of selec¬ 
tion is almost certainly more complicated in detail, the 
absence of additional information about the relation be¬ 
tween galaxy morphology and host halo properties makes 


ours a practical approach from a phenomenological per¬ 
spective. 

Impact of cosmology: the c-125 and c-250 simulations 
from which we drew our haloes assume h = O.T, Qm = 
0.29, Qb = 0.04T, as = 0.82, and Us = 0.96. These val¬ 
ues are slightly different to the best-htting Planck val¬ 
ues {h = 0.6T, Qm = 0.32, Qb = 0-049, as = 0.83 
and Us = 0.96; Planck Collaboration 2014), which 


themselves have some uncertainty. The underlying cos¬ 
mology (most signihcantly Qm, as and Ug) affects the 
log(concentration)-log(Mvir/M0) (c-m) relation of the 
haloes in the z = 0 snapshot, and hence the Vso values 
of their galaxies at hxed stellar mass. To investigate the 
robustness of our results to reasonable variations in cos¬ 
mological parameter values, we use the model of|D iemer| 
& Kravtsov (2015) to generate approximate c-m relations 
over the mass range of interest. In particular, we begin 
by htting a power law to the c-m relation generated by 
the cosmology of our simulations [pivoted at the approxi¬ 
mate median log(Mvir/M q) of the POT sample, Il.T], and 
then calculate the maximum fractional variation in the 
normalization when cosmological parameters are varied 
within 2(7 Planck uncertainties. We hnd that the nor¬ 
malization can be changed by at most ~ 10 per cent in 
this wayp^ To determine the effect this would have on 
our parameter constraints, we repeat the TFR MCMC 
analysis (for the Upeak model) with all halo concentra¬ 
tions shifted downwards by 10 per cent. We hnd that the 
best htting value of d is raised by ~ 2cr to — 

(c.f. Tablehrst row), indicating that a small amount 
of halo contraction is now permitted. The upper limit on 
the AM scatter is raised slightly to 0.26 (Icr). All other 
constraints are affected at the < 5 per cent level. 

Halo concentrations: estimates of concentration made by 
halo hnders may be unreliable, especially for haloes com¬ 
posed of few particles. We have checked that median con¬ 
centrations of ROCKSTAR haloes in the c-125 and c-250 
boxes agree to better than 5 per cent with the model 


of Diemer & Kravtsov (2015), and that the median frac¬ 


tional deviation between different ways of calculating 
the concentration is < 5 per cent (see also Section . 
We therefore conclude that there is at most a ~ 5 per 
cent systematic error in average concentrations, and have 
shown above (in the context of the impact of cosmology) 
that this has relatively little effect on our results. 

Stellar masses: we have taken stellar masses for the 
POT galaxies from the NSA, and the SMF for AM 


from Bernardi et al. (2013). Several prescriptions exist for 


The change in the slope of the c-m relation as cosmological 
parameters are varied within 2a Planck uncertainties causes 
fractional changes in the concentrations over the mass range 
of interest that are around 3 times smaller than the normal¬ 
ization difference, so may be safely ignored. 
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Figure Al. Comparison of the POT data (red points) to 1000 
mock TFR data sets generated using the maximum-likelihood 
parameter values of the Vpeak model (top panel; see also Ta¬ 
ble]^, and with parameters individually perturbed to illus¬ 
trate their impact on the TFR (other panels). The blue points 
are at the median velocities of the mock galaxies at the stel¬ 
lar masses of the POT galaxies, and the error bars show the 
interquartile range. 
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Figure A2. Same as Fig. |A1| but for the MSR. 



















































































































































































































































































































26 H. Desmond and R. Wechsler 


mapping observed luminosity to total luminosity to stel¬ 
lar mass, and their relative merits have been discussed 


extensively in the literature (e.g. Bell & de Jong 2001 


McGaugh||2005b Bernard! et al.||2013 ). However, a frac¬ 

tional shift in stellar mass that is the same for all galaxies 
of a given mass would not affect our results because the 
real and mock data would be affected in exactly the same 
way. Our constraints would only be modihed if the masses 
of the POT galaxies were affected differently to the aver¬ 
age for the entire galaxy population, which sets the shift 
in the SMF. Since the POT galaxies are typical in terms 
of colour and other properties (POT), this is unlikely. 


Impact of disc formation on halo: we have modelled the 
effect of disc formation on halo density prohle using the 
generalized adiabatic contraction model of [Dutton et al.| 
( 200T| . Although this model does allow for a halo ex¬ 
pansion (and should therefore capture the leading order 
effect of associated processes), the true mass dependence 
of the expansion may not be included. In particular, it is 
expected from some cosmological simulations that the ef¬ 
fect of stellar feedback on the dark matter halo is largest 
for log(M*/M 0 ) 8.5 ( |Di Cintio et al.||2014a|b ), whilst 
in our model the effect is largest at the peak of the stellar- 
to-halo mass fraction [log(M^/M 0 ) 10]. Changing the 

stellar mass dependence of the effect of baryons on the 
halo density prohle would alter the shape of the TFR (for 
example, the jPi Cintio et al.||2014b| model would reduce 
the predicted velocity more at the low end of the POT 
mass range, steepening the predicted TFR), but would 
affect none of our main conclusions. We note, however, 
that haloes produced in the EAGLE simulations (which 
also include the effect of baryons) are most contracted 
around the peak of the stellar mass fraction, an effect 
captured to first order by our generalized adiabatic con¬ 


traction prescription when u > 0 (Schaller et al. 2015 
Schaye et al.||2015 ). 


Mass-to-light ratios: we have assumed that the mass-to- 
light ratios of the bulge and disc are equal. This assump¬ 
tion is required to map the bulge-to-disc i-band fiux ra¬ 
tios observed in the POT study to mass ratios required 
for our rotation curve modelling. It is, however, generally 
accepted that galactic bulges have higher mass-to-light 
ratios than discs; thus the bulge should be more massive 
relative to the disc for a given fiux ratio. We have in¬ 
vestigated the effect on Vso of variations in the ratio of 
the mass-to-light ratios of the bulge and disc from 0.3 to 
3, and find that changes are typically at or below the 1 
per cent level, due to the relatively small bulges in these 
late-type systems. 


Velocity dispersion: we assume that the model galaxies 
in the population from which the POT sample could have 
been drawn have negligible velocity dispersion relative to 
their rotation velocity. It is expected that galaxies with 
properties similar to those in the POT sample have ve¬ 
locity dispersions of order 20kms“^ (McCaugh, private 


communication), which would cause a relatively small 
shift of our model galaxies towards lower Vso values. In 
addition, we assume our model galaxies are thin when we 
simulate their rotation curve; real discs have finite thick¬ 
ness which reduces their rotation velocity at fixed mass 
and size. For realistic thicknesses (e.g. h/Ra ~ 0 . 2 ), this 
is a < 5 per cent effect (e.g. Casertano||1983| ). 


Dynamical effect of gas: we assume that the cold gas 
in galaxies contributes negligibly to the total dynamical 
mass, and hence to Vso- Using the approximate scaling 
relations for the gas and stellar masses and disc scale- 


lengths presented in Dutton et al. (2011), we expect the 


gas disc to increase typical galaxies’ Vso values by at most 
around 2 per cent in the mass range we are considering. 


Uncertainty in stellar masses: we have assumed no un¬ 
certainty in any stellar mass estimate. If this uncertainty 
were non-zero, each AM scatter value that we quote 
would become the quadrature sum of this uncertainty 
and the intrinsic scatter associated with AM. Since it 
is the true rather than observed stellar mass which de¬ 
termines a model galaxy’s rotation curve, however, one 
would need to perturb the stellar masses output by AM 
according to their uncertainty before continuing with the 
calculation of rotation velocity and disc scalelength. This 
would increase the intrinsic scatter of the theoretical TFR 
and MSR. However, the effect would be small because 
the low stellar-to-halo mass fractions induced by AM (at 
most ~ 5 per cent) make a galaxy’s own contribution to 
its rotation curve far less than that of its halo. The slope 
and normalization of the relations would be on average 
unaffected, as would the correlation of the TFR and MSR 
residuals, which is recorded at fixed stellar mass. 


Halo sphericity: we assume haloes are perfectly spher¬ 
ical. However, cosmological simulations create triaxial 
haloes with smallest-to-largest axis ratios of the order 


of 0.5 (Allgood et al. 2006 Despali, Ciocoli & Tormen 


2014). This is expected to be increased by the process of 
galaxy formation ( Kazantzidis, Abadi fc Navarro|| 2010 ), 
and weak leasing observations provide only weak evi¬ 


dence that haloes are not spherical (van Uitert et al. 


2012). Triaxiality would cause the rotation velocity of the 


disc to be a function of angle in addition to radius, intro¬ 
ducing some thickness to the observed rotation curve and 
increasing the scatter in the TFR, but the magnitude of 
the effect is unlikely sufficient to affect our conclusions. 


Halo density profile: we assume that haloes have a density 
profile perfectly described by the NFW form. Deviations 
from pure NFW would cause variations in the circular ve¬ 
locities of the haloes for given rs,kiypin and rvir, generating 
additional scatter in the TFR. However, this approxima¬ 
tion is unlikely poor enough to constitute a significant 
source of systematic error in our analysis. 


Disc intrinsic axis ratio: we take all discs to have an 
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intrinsic axis ratio of 0.2, as assumed by POT. It is pos¬ 
sible, however, that these could be signihcantly higher, 
especially at low mass (Geha, private communication). 
Underestimating the axis ratios would cause the inclina¬ 
tions to be underestimated, and hence the inclination- 
corrected Vso to be overestimated. To explore the impact 
this could have on our parameter estimation, we repeat 
the MCMC analysis of Section with two alternative as¬ 
sumptions. In the hrst, the intrinsic axis ratios of all the 
POT galaxies are set to 0.6. In the second, we make the 
axis ratio a linear function of log(M*/M0), falling from 
0.6 at log(M*/M0) = 9 to 0.2 at log(M*/M0) = 11.2. 
(This modifies the slope of the observed TFR as well as 
its normalization.) In both cases, we hnd that our param¬ 
eter constraints are affected at the < 5 per cent level. 


• Rotation curve sampled at 500 logarithmically uni¬ 
form intervals over 10“^ < r/Rd < 30. 

• Number of mock data sets for likelihood evaluations: 

200 . 

• AM scatter resolution: 0.02 dex. 

• v resolution: 0.1. 


APPENDIX C: TECHNICAL DETAILS OF 
THE METHOD 


In order for the likelihood evaluations to be sufficiently 
rapid for a full MCMC analysis to be feasible, the halo 
AM step cannot be performed each iteration. We there¬ 
fore began by creating and storing arrays of the M* val¬ 
ues associated with each halo for each AM parameter 
and scatter, in addition to the haloes’ Mvir, concentra¬ 
tion and A values, and (for Section 5.4) the ROCKSTAR 
‘pid’ variable that distinguishes satellites from centrals. 
We then created a look-up table to map from a point in 
{Mvir, concentration, M^, A} space to a point in {Rd, 
Vso} space using the method described in Section]^ par¬ 
titioning each space into a grid with the properties de¬ 
scribed below. This was done separately for each value of 
the v parameter, which affects this map. Then, in the like¬ 
lihood function and for a given set of parameter values, 
we determined the number of haloes in each {Mvir, con¬ 
centration, M*, A} grid cell, and incremented the counts 
in the TFR {M*, Vso} or MSR {M^, Rd} arrays by the 
corresponding amount. These arrays were then used to 
generate mock data sets. This procedure allowed each 
likelihood evaluation to take only 25 s. More quantita¬ 
tive information concerning the convergence conditions 
and resolution of our algorithm is given below. 


• Convergence requirement for iterative calculation of 
disc scalelength (see Section]^: 2.5 per cent. 

• Convergence requirement of generalized adiabatic 
contraction: 5 per cent. 

• Resolution of grid of log(Mvir/M0), 
log (concentration), log(M*/M0), log(A) values: 50 
cells in ranges 9.5 : 13.0, 0.4 : 2.1, 9.0 : 11.2, —2.3 : —0.4 
respectively. These limits were chosen to enclose > 99 
per cent of model galaxies for all AM parameters. 

• Resolution of log(V8o/km s“^) grid: 50 cells in range 
l.TO : 2.T5. 

• Resolution of log(Rd/kpc) grid: 50 cells in range 
- 0.8 : 2 . 2 . 
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